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A compact gear box is one of the main considera- 
tion in the economics of a machine tool design. Compactness 
is deuendent upon the layout and speed diagrams. A mathema- 
tical model for the compactness of a 2 - speed gear box is 
formulated. Interior Penalty Function approach is used to 
arrive at an optimum speed diagram and the gears. The re- 
maining components are designed and inspection data for gears 
are generated for the shop floor use. 
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CHAPTER I 


INTRODUCTION 

The designer of a multipurpose machine tool has to 
provide a certain speed range which covers the requirements 
of the different operations, types and shapes of work-pieces 
and qualities of the surfaces that are to he machined. These 
requirements together with the material of work-pieces and 
tools determine the most economic speeds for the two machi- 
ning movements: the cutting movement and the feed movement. 

The ideal situation for handling the two machining 
movements would he to provide stepless spindle and feed 
speeds. However, the majority of machine tools rely on she 
pped drive mainly because of rigidity considerations and 
ease in manufacturing them. Thus, the gear box, 1 giving 
the required spindle and feed speeds is one of the most 
important sub-assembllos . 

There are several possible arrangements by which the 
required number and magnitude of spindle speeds may be achie- 
ved for the same speed of the input shaft. These arrangements 
are called the ray, tree, or layout diagrams (Pig. 1.2'^ and 
Fig. 1.4). For each such layout diagram there are several 
ways by which the desired spindle speeds may be obtained. 

This is possible because the transmission ratios between the 

T7 Fig. 1.1 £nd Fig. 1.3 show the schematic arrangement of 
shafts and gears in developed view for possible 10 and 
18 speed gear boxes. 
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different shafts nay "be interchanged or altered* Any such 
arrangement of obtaining the required spindle speeds is 
depicted in the speed diagram. Therefore the problem of a 

gear box design is complete once a suitable ray and corres- 

1 

ponding speed diagram is chosen and the components are 
designed for the adequate power transmission at the speeds 
obtained by the speed diagram. 

1*1 Critical Review of Past Investigations! 

F» Koenigsberger^ 1 has put some criteria for the 
suitable choice of a ray diagram. But the author is not very 
decisive about the appropriate choice of speed diagram for 
the proposed ray diagram. The dynamic performances of a gear 
drive, according to him, depend largely upon the bending and 
the torsional deformation of shafts, gears etc. The stiff- 
ness of a transmission between input and output shaft decreases 
with increasing lengths of the shaft and the ‘Centre distance 
between the shafts. Therefore, the author agrees in principle 
that such speed diagrams should be chosen which will require 
least centre distance. 

G. White proposes a method to select the gears 
for a four speed transmission with two composite gears. The 
arrangement for which maximum number of similar gears can be 

1. Henceforth mention about speed diagram shall always 
be corresponding to the ray diagram. 
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used is also discussed by the same author. The algorithm is 
to solve a set of non linear algebraic equations with number 
of teeth in different gears as the variables. This works well 
with low number of transmissions from a gear box but the equa- 
tions become very complicated and non linear in nature with 
the increase in the transmissions. 

In a machine tool manufacturing concern, the gear box 
is designed by arriving at a workable ray and a corresponding 
speed diagram based upon the designer’s past experience. The 
general norm which is followed in the industry states - 
’’Out of innumerable speed diagrams available, those should be 
chosen for which the left side of the speed diagram is most 
concave”. It follows from here that the intermediate shafts 
should have as high speeds as possible. 

■ 11 ) 

F. A. Cinadr has recently solved a gear box optimi- 
zation problem which is an eye-opener for the engineers using 
computers in machine design. With the help of his programme 
one can find the optimal dimensions of machine elements like 
gears, clutches, bearings, keys etc. for a most compact 12 
speed gear box. But the programme suffers from the following 
disadvantages; 

1. The modules have not been selected as independent 
design variables. 
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2. The programme can handle the design of a l2_ speed 
gear box for one particular layout diagram only. 

The main objection to the arbitrary selection of a 
speed diagram is that ono may not choose the speed diagram 
which gives the minimum centre distance between the shafts 
of the gear box. Therefore, in the view of the present 
author, if the speed diagram is arbitrarily selected then the 
optimal choice of the remaining gear box components has little 
significance . 

1,2 Present Investigation - Its Scope: 

Here a mathematical model of a z-speed gear box has 
been developed. Depending upon a selected optimality crite- 
rion, corresponding to any layout diagram, an optimal speed 
diagram can be obtained. All the members of the gear box are 
designed for all loading conditions out of which the best 
dimensions are chosen. For this, a computer programme, to 
avoid excessive hand computations has been developed* Impor- 
tant inspection data for checking these components is also 
generated in the same computer programme. 

In view of the fact that the compactness of the gear 
box leads to economy of resources in terms of saving of 
material, machining time and weight, the minimum centre dis- 
tance between the first shaft and the spindle is chosen as 
the optimality criterion. 



5 


The optimal design problem is formulated as a mathe- 
matical programming problem. The optimization algorithm 
chosen is the Fletcher and Powell's variable metric method. 
The method is iterative in nature and seeks the optimum from 
within the constraint set. The method is discussed in detail 
in Chapter IV. 

1.3 Definitions? 

Gear Box: A sub-assembly of the machine tool con- 

sisting of spindle, motor and other 
intermediate shafts with gears mounted on 
them. Through different pair of gear connections the desired 
output speed of the spindle is obtained. (Fig; 1.1 gives the 
schematic diagram of an IS speed gear' box). 


Layout Diagram: This indicates the mode of transmission 

(tree or Ray 

diagram) of power and speed from motor shaft to 


spindle shaft for a given number of 
spindle speeds. (Fig. 1.2 represents a layout diagram for 
an 18 speed gear box). The lines (U . U-^ etc.) indicate 
the transmission from one shaft to the other. The distance 


between two such transmissions is called the step of trans- 

1 S 

mission ratio (As 0 ± on shaft No. 2 or 0 on shaft No. 5). 


Group of If the spindle speeds a of a gear box 

transmission: 

be split in K simple numbers (viz. for 
z= 18 =3x3x2, K=3 for z = 24 
= 2x2x2x3or4x3x2,K = 4or3 respectively 
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depending on how the designer wishes to split them) then there 
are K groups of transmission. (For the arrangement shown in 
Fig. 1.2 there are 4 groups viz. 18=3x3x1x2). 

Step of Spindle This indicates the common ratio of the 
speed; 

series of spindle speed in case of geome- 
trically stepped speeds. This shall he indicated by "0". 

Group This is the power index of "0" of the 

characteristic ; 

common ratio of speeds from a group (viz. 
in Fig. 1.2 s l for group 1, 9 for group 4). 

Speed Diagram; This represents the actual speed for all 

the shafts together with values of diffe- 
rent transmission ratios for a particular layout diagram. 

The shafts are plotted as equidistant parallel lines and the 
speeds on a shaft are plotted on logarithmic scale. (Fig.6.l) 

Speed range This is the ratio of maximum to minimum 

ratio (B); 

spindle speed. 

Composite gears; When a driven gear acts as a driver gear 

for the next shaft, it is called a com- 
posite gear. 

Back gear shafts If nth shaft in a drive be the back gear 

shaft then it enables (n + l)th sha.ft 

<n 

to have all the speeds of the (n + l)th shaft with reduc- 
tions which is not possible by a single pair of gears. 



Main drive of 18 speed gear box 



Asterisked gear and their module indicate 
desigin variables. 


Fig. I -1 



o 








CHAPTER II 


DESIGN PHILOSOPHY 

In this chapter, important factors which aid in the 
choice of the number of spindle speeds and corresponding 
selection of layout and speed diagrams are enumerated. The 
total number of gears necessary to obtain required spindle 
speeds is not unique. Some qualitative considerations are 
presented for the selection of the number of gears. 

2.1 Step Ratio j 

The ratio of two consecutive speeds of spindle is 
called the step ratio. Most economic cutting speed for a 
particular material and job diameter of a work-piece is unique. 
In order to machine a job of any diameter and material most 
economically, stepless spindle speed is necessary. But from 
the considerations of rigidity, economy and ease of manufac- 
ture, the stepped drive is always preferred to the stepless 
one . 


The spindle speeds may be arithmetically, geometri- 
cally or logarithmically stepped. With the stepped' drive it 
will be necessary to establish upper and lower limits for the 
economic cutting speeds. A geometric progression for spindle 
speeds is accepted because it can handle a large diameter 
range for the same speed range and is superior from the design 
and kinematic aspects. Since the intermediate transmissions 
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are also stepped in geometric progression, it has the added 
advantage of providing easy calculations. The number of spin- 
dle speeds is kept as low as the economic handling of large 
numbers of diameters and materials tllow. This requires a 
finely stepped feed drive which is more economical as compared 
to fjnely stepped spindle speeds. For a z - speed gear box, 
the common ratio 0 is obtained for a given speed range ratio 
E by the relation i 
z-1 

0 = B (2.1) 

where 

t _ Maximum spindle speed 

b ~ Minimum spindle speed 


Rewriting Eq. (2.i) in a different form 


z-1 


or 


log B 
log 0 


= 1 + iogB 

log 0 


( 2 . 2 ) 


Maximum and minimum spindle speeds and step ratio are not 
chosen arbitrarily. German Standards DIR 804 and DIN 805 
provide standards for both spindle and feed speeds. These 
Standards are based on the fact that they cover the majority 
of economic cutting speeds depending upon the number of spin- 
dle speeds . The plot between B and z with constant 0 values 
is shown in Fig. 2.1. From the graph the designer can weigh 
up the various solutions at his disposal for economic selec- 
tion of number of spindle speeds. 



is 


2.2 Layout Diagram} 

If the number of spindle speeds z is split into K 
simple numbers then there are K groups of transmissions in 
the gear box. Bach simple number represents a group and its 
magnitude indicates the number of transmissions from that 
group. For example, for a 24 speed gear box 

z = 24 = 1 04 . .2 ft] . 3 . 4 iff] (2.3) 

Here there are four groups in the gear box with number of 
transmissions from each group as 2, 2, 2 and 3 respectively. 

The group arrangement shown in expression (2.3) conforms to. 
standard practice. 

Splitting of the spindle speed z may not always con- 
tinue to the last prime number. It depends on the number of 
transmissions that the designer wants to allow in a group. As 
a matter of fact in obtaining the number of groups, if the 
factor H 2 n appears more than once, it should be combined to 
give four transmissions from one group instead of two trans- 
missions from two groups. For example, for a 24-speed gear 
box the arrangement shorn in expression (2.3) should be avoided. 
The obvious reason is the economy of space by saving a group. 
However, according to standard practice, the maximum number 
of transmissions from a group seldom exceeds 4 because in this 
case the shaft length increases with decreasing rigidity. 
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The total number of ways in which a layout diagram 
can he chosen is limited. If out of K groups D have equal 
number of transmissions, then the total number of layout dia- 
gram M is given by 

M «, ' K l Kl (2.4) 

A d : 

For an 18 speed gear box, one of the possible arrangements is 
z = 18 = 1 [sj . 2 [sj . 3 [sj 
Here the number of groups K is equal to 3, and number of 
groups D, having equal number of transmissions (1st and 2nd 
groups) is 2. Therefore, the total number of possible layout 
diagrams is 

M = — H - L = M 

Some of these 18 layout diagrams are shown in Fig, 2.2. 

The above discussions about total number of layout 
diagrams are valid only when there is no back gear shaft at 
any stage of the gear box. Fig. 3.3 shows a back gear connec- 
tion between the shafts of a gear box. The advantages of 
using back gears are twofold. Firstly, the transmission 
ratios which cannot be obtained by a single gear^pair can 
be obtained here. Secondly, if the nth shaft be the back 
feear shaft then the arrangement enables (n + l)th shaft to 
have all the speeds of (n - l)th shaft. This is done simply 
by direct coupling. 
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To incorporate the consistency in the group nomen- 
clature for computer programme, two groups are considered' 
whenever a back gear is present, as shown in Fig. 3.2, namely 
the groups containing the transmissions from (n - l)th shaft 
to nth. shaft (back gear shaft) and from nth. shaft to (n + l)th 
shafts 

It is only at the later stages (towards the spindle) 
that a back gear arrangement can be used with advantage be- 
cause then they provide for increased number of speeds. In 
the. presence of back gear shafts, total number of possible 
layout diagrams is difficult to determine. 

The arrangements which have more than two trans- 
missions in the last group(i.e. on the spindle) should be 
rejected because of excessive dynamic loading on the spindle. 
For example, for a 24 - speed gear box, arrangement, such as 
24 = 1 [ 4 ] , 2 [ 2 ] , 3 should be rejected. Number of 

transmissions from a group sued their occurrence in the group 
order are factors which account for the size of the gear 
box. Group orders which have larger number of transmissions 
in the lower group (towards the 1st shaft) than in the 
higher are preferred. The gears in the lower groups are 
nearer to the motor shaft as compared to the gears in the 
higher group. Conf >quently, gears in the lower group run at 
comparatively higher speeds and transmit less torque. There- 
fore, they are smaller in size. So for compactness of the 



gear box, more number of transmissions are incorporated in 
the lower groups- In certain exceptional cases this rule 
is violated. For example, when a designer has to select 
a layout diagram for 12 speed gear box, incorporating flexi*r 
bility in the design to convert it to 18 speed gear box, with 
minimal change in the assembly, the following arrangement is. 
preferred: 

12 = 1 jVj . 2 \f] . 2 [2] 

In this by adding one more transmission to the first group 
It is converted to an 18 - speed gear box with following 
arrangement : 

18 = i [ 3 ] . 2 [3] . 2 [2J 

Layout diagrams with smaller variation of speeds on 
the intermediate shafts are preferred. With such selections, 
a pair of gears in a group, which is optimally designed for 
lowest speed does not become much over designed for higher 
speeds of same shaft. It follows that if x , Xg, ... x^ be 
group characteristics from 1st to the kth group, they should 
satisfy the following inequality : 

x l< x 2<“*< x k C2 * 5) 

In some cases this inequality is violated. One such example 
is the use of two composite gears in the second group. In 
Appendix - 2 the arrangement where x 1 > Xg is superior to 
the one where x / x is shown. 

X ^ (-* 



2.3 Number of gear pairs required: 


17 


At this stage, the total number of gears necessary 
to meet the speed requirement has to be decided upon. The 
number of gears in the transmission can be reduced by using 
composite gears, where one driven gear in a group acts like 
a driver gear for the next group. This arrangement results 
in increased centre distance but the number of gears is re- 
duced. 


A pair of gears can be connected either by clutches 
or by sliding the cluster gear on the splined shaft. The 
former is the constant - mesh - type arrangement which is 
more silent and can use helical gears but is expensive. The 
sliding block or cluster gear arrangement is used at high 
torque and is cheaper but noisier. 

2,4 Inter - departmental activity: 

The selection of intermediate shaft speeds is not 
unique. The selection has to be a compromise between various 
conflicting situations. Some cf those are stated below. 

(a) Speeds of the irtermediate shafts should be 
as high as possible to reduce the torque on 
the shaft. Transmission ratios less than 0.25 
and greater than 2.0 are seldom allowed. 
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(1) With high sp^ed of gears the dynamic loading 
due to inaccuracy in the machine cut gears 
increases tremendously and, therefore, high 
quality gears have to he used. Pitch line 
velocity should not exceed 25 m/sec. in majority 
of machine tool gear drives. 

(c) Bearing loading conditions favour lightly loa- 
ded high speed shafts than heavily loaded low 
speed ones. But then, if the speed is more, 
the gears are smaller in size and the load on 
the shaft increases. 

(d) Final speeds of the spindle must maintain a 

2 % accuracy from their respective desired speeds. 

It is evident that the sizes of components like gears, 
shafts, clutches depend upon the choice of the speed diagram. 

So the selection of a speed diagram is the vital decision an 
engineer has to take before proceeding further. 

The choice of the best speed diagram depends upon 
any desired criterion of optimality. For most practical pur- 
poses economy should be the guide line. To get the most 
economical gear box all the factors that influence economic 
production, such as the dimensions of all members, their 
respective materials, the manufacturing time, the processes, 
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man. and machine hours involved have to he considered. This 
becomes a formidable task. A simpler optimality criterion 
such as compactness greatly reduces the task of quantifying 
some of the above influencing elements. 




FIG. P.i 




Ztr 




CHAPTER III 


THE MATHEMATICAL FORMULATION 

In this chapter, the design variables used in the 
mathematical modelling of the gear box are introduced and the 
objective function is formulated. The design constraints 
are expressed in the form of inequalities and a penalty func- 
tion to account for the constraints is generated. 

3.1 Design Variables: 

Once the number of so indie speeds and that of groups 
are decided upon, an appropriate layout diagram remains to 
be chosen. The layout diagram indicates the number of speeds 
and the ratio between any two speeds of a shaft. Therefore, 
if the smallest (or any other) speed of a shaft is known all 
the other speeds of that shaft can be calculated. Given the 
input speed and the spindle speeds, lowest speed of interme- 
diate shafts can be arrived at in such a fashion that lowest 
transmission ratios from all the groups meet the spindle 
speed requirements. These transmission ratios can be changed 
to give infinite number of altogether different speed diagrams 
for arriving at specified spindle speeds, thereby indicating 
that some of these are the independent design parameters. 

As long as. the teeth ratio is equal to the trans- 
mission ratio between two shafts, any number of teeth and 
module value for pinion and gear can be used for transmission 
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in that group containing the shafts. However, the gears 
must he of adequate strength to transmit the power to the 
spindle effectively. The transmission ratio between a pair 
of gears in a group, the number of teeth of the pinion and 
its module determine the centre distance between two shafts 
of that group. Therefore, whenever these quantities can be 
changed arbitrarily they form a set of independent design 
variables. If in some groups, the number of teeth of one of 
the pinions cannot be selected independently and the trans- 
mission ratio between two gears is not a free choice, then 
there is no independent design Parameter in that group. 

3.2 Nomenclature of Important Components; 


The gears are named in double subscript, Z^ . 

First subscript indicates the number of the shaft on which 
the gear is mounted and the second subscript indicates the 
number of the gear. A definite order is followed in the se- 
quence of numbering the gears. The driver gears are numbered 
first in the order of increasing transmission ratios. Then 


the driven gears follow a similar order. Corresponding to 

every gear there is a notation for transmission ratio TJ , 

modules m. minimum speed S.., maximum speed T^, and com- 

th ^ th ^ 

posite gears L..» The 3 speed of the i shaft is denoted, 
by R. .. Some of the number of teeth, modules and transmi- 

J 13 

ssion ratios are arbitrary choices (in the bounds of the design 


constraints) and are represented by a vector X of design 
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variables. The remaining are determined as their non - 
linear combinations. 

3.3 Design Space: 

If the number of independent variables in a parti- 
cular configuration of a gear box is n, then the n-dimensional 
Eucledian space is known as the design space. Every point 
in this apace represents a design even if it is absurd (as 
negative number of teeth or a gear digging into the other 
shaft). The movement from point to point in this space re- 
presents a change in design. Any point in the space is repre- 
sented by a vector 

x = U, V ... , xl ‘ (3.1) 

3.4 Constraint Surface: 

All the points in the design space do not represent 
a feasible design. Each of the design requirements will be 
satisfied in some region of space and violated in the other. 
The surface that separates these two regions is called a 
constraint surface. Every design requirement has a corres- 
ponding constraint surface. The region of space in which all 
design requirements are satisfied is called the feasible 
region or the acceptable region. 

Each of the constraints is represented by 
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g ± ( X ) ^0 , i = 1, 2, ... N c (3.2) 

where 

N c = number of constraints. 

3.5 Objective Function 

The function f ( X ) of design variables, the value 
of which is a measure of the merit of the design and which 
must be optimised in the feasible region is called the ob- 
jective fu'iction. The best feasible design is one for which 

g. ( X ) ) 0 , i = 1, 2, ... N c (3.3) 

. and 

f ( X ) is optimum. 

In the present work the sum of the individual center distan- 
ces is chosen as the optimality criterion. 

The gear boxes can be divided Into four general 
categories depending on their complexities. Gear boxes may 
have 

1. No composite gear in any group. 

2. One composite gear in a group 

3. Two composite gears In a group 

4. Back Gear at some stage. 

There could be other classes as well, which can be 
deduced from the above categories. 
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The following factors are common to all types of 


gear boxes in objective function formulation. 

1. The set of design variables correspond to some of 

the elements in matrices fz^"} , [ m ij3 5 ^ijl * Tliere f ore 
any change in design variables is accompanied by a change 
in the corresponding elements in the matrices. The other 
elements of the matrices are combinations of design variables, 
Hence, they get modified every time the design vector is 
changed . 


2. The objective function in all the cases is 


whe ri 


f ( X ) 


1=1 


_.th 


m = number of groups, 

. = Center distance in the i v '“ group 

Z il * ( U il + 1 ‘ 0) * m il 


2 . U 


il 


(3.4) 


th 

3. Once the lowest transmission ratio from the i group is 
th 

known, the j transmission ratio can be calculated as 


U 




= n u • 0 


(j-l).x i 


(3.5) 


isli 

x^ is the group' characteristic of the i group.. 


whe re 
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4. Fvery speed of the (i-l)th group generates as many as 

n i-l s P efids on shaft, where n^_^ is the number of 

transmissions from the. (i-l)th group. If represents 

the number of speeds of (i-l)th shaft, then total number of 
th 

speeds of i shaft is 


The input shaft has only one speed R . The speeds of the 
second shaft are obtained by the relation 

Rg ^ = * ^lj * 5 = l}^j ... n^ (3.6) 

For every speed of second shaft, the speeds of the third 
shaft are given by 


R 3jk *2i * U lk ! 


., k = 1, 2, . . . ng (3.7) 

Evidently the minimum speed of the i^* 1 shaft is given by 
the relation 

R il = R 11 * U 11 * U 2l *** U i-l,l 


5. The number of teeth of all the gears in the i^ R group is 
determined if and the transmission ratios , 5 = 1, 

2, . . . n^ are known. This is given as 




(l.o + u i3 _) . u i3 
TT.o +'u^) 


2 =1, 2, . . *n i 

(3.8) 


6. In every case the lowest transmission ratio from the 
last group is not arbitrary. If m be the total number of 
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groups, then it follows from the equation 3.7 


Vi,i - hi ' °n 



u 


m-1,1 



Therefore , 


R 


U 


m+l,i 


ml 


(3.9) 


E H • U 11 * U 2i .*•* u m-l,l 


where 


R m+l l ^ s ' t ^ e l° wes_t given spindle speed. 


It follows that C i t s can he determined If 
and are known for the i^ h group. Methods of finding these 
quantities for every group changes from one configuration of 
the gear box to the other. Actually depending upon the con- 
figuration of the gear box, the number of independent design 
variables changes. This is discussed in the subsequent 
paragraphs for the earlier mentioned types. 


Type 1; A four speed gear box with eight gears is shown 
in Fig. 3.la. None of the gears in both the 
groups is a composite gear. For a given mini- 
mum speed R q and the given input speed R , 
the minimum speed Rg^ of shaft 2 can assume any 
value by varying U arbitrarily and setting 


R 


U, 


21 


31 


R . U ni 
11 11 


(3.10) 
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The gears Z and Z and their modules can 
11 2l 

also he arbitrarily chosen without affecting 
the compatibility. Therefore, 'the number of 
independent design variables in this case is 5. 

Type 2 : The same 4 spindle speeds have been achieved by 
7 gears in Fig. 3. lb by using one composite 
gear in the second group. In this case also, 

U can be arbitrarily varied and is found 
from the Fq. 3.10. The number of teeth and 
module of Z^ in the first group can be chosen 
arbitrarily but the number of teeth and module 
in second group is not a free choice and is 
automatically determined from group 1. There- 
fore, the number of independent design variables 
reduces to 3 in this case. 

Type 3 : Fig. 3.1c represents a gear box arrangement to 

reach the same four speeds of the spindle 

using only 6 gears. Both the driven gears of 

the first group are composite gears. It has 

been shown that once the number of teeth and 

module of Z_^ are arbitrarily selected, all 

the numbers of teeth can be determined as a 

function of 0 , R„ and R__. This means that 

ol 11 

there can be only one possible value of trans- 
mission ratio for given Therefore, 
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Type 4s 


the number of independent design variables 
reduces to 2. 


The expression for U in this case is 

1 - H . 0 . (0 + l) 

U = ■ 

11 0 


(3.11) 


Fig. 3.2a indicates the 4 speeds on the 4 
shafts by a back gear arrangement. The centre 
. distance between 2nd and 3rd is the same as 
centre distance between 3rd and 5th shaft. 

The arrangement uses 8 gears and its layout 
diagram is shown in Fig. 3.2b. The transmi- 
ssion U is not arbitrary because the speed 
it generates on the second shaft is trans- 
ferred to the 4th shaft directly by clutch and 
has to be one of the desired speeds. Any one 
of the transmission ratios and is 
arbitrary. Others are given by the relations 


R. 


U 


43 


11 


‘11 


R 


U, 


41 


31 


(3.12) 


B 11 • *U 


U 


21 


The number of the teeth and module of Z 11 and 

Z_ are arbitrary. Z can be written as 

ol 
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z 3l (1.0 + U 3l ) : mg£ 
whe re 



(3.13) 


The number of design variables is, therefore, 
5 in this case. 


The discussion of the above four types has been 
included for the sake of completeness. The computer program 
by itself finds the number of design variables and evaluates 
the objective function for any of the above types of their 
combinations . 

3.6 Normalised Constraints'. 


The general form of expression of design constraints 
is 

g ( X ) = 0.01 - f x ( X ) >0 

X 

g 2 ( X ) = f 2 ( X ) - 10,000^0 

where it is required that 

f x ( X ) ^ °- 01 

* f 2 ( X ) 10,000 

The above is an example where the magnitude and 
sensitivity of the two constraints with respect to change in 
design variable ate quite different. This results in 



32 


difficulty to arrive at the optimum point efficiently. The 
disparity in magnitudes and sensitivities of the design cons- 
traints can he avoided by normalising them, i.e. by forcing 
them to take values between 0 and 1. If the inequality cons- 
traint expresses the difference between two variables, then 
it should be normalised with respect to the variable of 
greater value. Even when the magnitudes of the constraints 
are controlled, their sensitivities may vary widely. Never- 
theless, normalising always Improves the handling of 
constraints. 

3.7 Equations of Constraint Surfaces; 

It has been stated that the number of independent 
design variables and the number of gears required to obtain 
the desired spindle speeds vary with selection of layout 
diagram and the arrangement in gear box. Therefore, the total 
number of constraints is problem-dependent. 

3.7.1 Positivity of Design Variables; 

Negative values for the number of teeth, transmi- 
ssion and the modules of the gears are meaningless. To 
restrain the minimisation algorithm from going into this 
region of design space, positivity constraints are added. 



S3 


( X ) = 



•\r 'N 

' i 


^max"'*"" 


A . ),^ 

min l 


\0 . i - 1>2 


where 

n - Number of independent variables, 
th 

X ± = i design variable, 

(X min ) = Minimum value of i th design variable, 
(X maz ) = "maximum value of i th design variable. 


» n 

(3.14) 


3.7.2 Gear T'^eth: 


The minimum number of teeth of a gear blank should 
not be less than 14 in order to avoid excessive undercutting. 
This constraint is applied to all the gears in the drive and 
is expressed as 


g ( X ) = 1.0 - if — 3 = 1, 2, ••• N g 

i = n+l, n+2,...n+N 

o 

(3.15) 


3.7.3 Pitch Line Velocity: 


The dynamic loading of a me chine -cut gear increases 


tremendously with the increase in pitch line velocity 
of the gear. The r.p.m. in calculation of is obtained 


from the speed diagram. 

Z . .. m. . T 

l l : 

c = 


20,000 


, i = 1, 2, ... N 


g 


(3.16) 


The corresponding constraint equations are 

G. 


R ( X ) = 1.0 - iL - 

g i - 25.0 


0 j— 1,2, • • • N 


i = n+N +l,...N+2N 

CP ' 


g (3.17) 


g 



3.7.4 Wear Strength; 


All the gears in the drive should he checked for 
"bending and wear at rated horse power of the gear box. In 
view of the fact that amongst different types of gear 'failures, 
failure due to wear is most common, it is sufficient to check 
for the wear strength of the gears only. The expression for 
the horse power capacity in wear of a gear teeth is given by 


where 


(H.P. ) = X 

wear c 



1*8 

b , m . S 


z 


IU3 

10 6 

(3.18) 


X^ = Speed Factor, 

Y z = Zone Factor, 

S c = Service Factor, 
b = Face width of gear 
S = Speed of rotation 
m = Module , and 
z - Humber of teeth. 


The speed factor X is a function of r.p.m., s and 

c 

the number of working hours per day, H. This has been re- 
presented in a graph for different speeds and working hours 
in I.S. 4460-67. Empirical relationship is established here 
through the use of minimum deviation technique, and values 

of X are obtained upto 1$- accuracy from the relation 

C O 218 
0 . 8442 ( H - 0.821) * -1.385J 

= S e (3.19) 


X 
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The zone factor Y is a function of number of 

z 

teeth of the driver gear, Z end of the driven gear, Z , 

P ’ g 

This has been obtained from I.S. 4460-67 by interpolation. 

2 

Y z = A + B Z g + C Z (3.20) 

whe re 

2 

A = -1.3 + 0.0294 Z„ + 0.000169 Z 

P P 

E = 0.222 - 0.00204 Z + 0.0000115 Z 2 

P P 

C = — 0 , 00461 + 0.0000653 Z p + 0.0000003 Z p 2 

The service factor S c is defined for different 
combinations of materials for gear and pinion. The value 
here is taken as 1.125 for 0.4% Carbon Steel. 


The series of constraints g^' s are written as 


s i 


) = 1.0 


H.P 

(H.P. wear 


•y^O, j - 1)2, . . 

3 i = n+2Ng+l 


N g (3.21) 
. . »n+3N 

g 


3.7.5 Minimum Gear Blank Size 


Diameter of bore of the gear should be atleast 
2.35 times less than the pitch diameter to provide for 
sufficient strength of the annular portion left in the gear. 
To account for keys, a further factor of safety of 1.25 is • 
required. Therefore, maximum allowable bore diameter D^ in 

4 “ 

the i gear blank is given by 



Z. . m. - 2.35 . m. 
11 

D i = " L25"” 

or 

D ± - 0.8 . m i ( Z j _ - 2.35 ) 


and the corresponding constraint equations are 


se 

(3.22) 


where 


g. (X) = 1.0 - ^0, 3 = 1, 2, ... N g (3.23) 

i = n+3N +l,...n*4N 
g 7 g 


(dia). is the shaft diameter holding the i gear. 

d 


3.7.6 Modules: 


The modules of gears in most generally used machine 
tools vary from a lower limit of 1.75 to an upper limit of 
10.0. This presents 2 x m number of constraints where m is 
the number of groups in the drive. 

g (X) = 1.0 - -^T^O, j = 1,2, ..., m (3.24) 

3 

i = n+4N+i . . .n+4N +m 

S g 

and 

m . 

g ± ( x ) = 1.0 - > 0,3 = 1 , 2 , ... m ( 3 . 25 ) 

i = n+4N g +m+l, , . .n+4Ng+2m 

3.7.7 Transmission Ratio: 

There are constraints on the minimum and maximum 


values of transmission ratio from a group. Low transmission 
ratios lead to excessively big gears. Designers recommend 
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the upper and lower values of transmission ratios as 2.0 and 
0.25 respectively. This also adds 2 x m more number of cons- 
traints to the problem. These are written as 


and 


where 


% 


= 1.0 


g i °° 


l.o - 


■ ^°5 j = 1» 2 ? ••* m 

u n 

i = n+4H +2m+l 

g 

(N j - 1) .x, 

U . 0 3 3 
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2.0 


>0 


(3.26) 
,n+4N +3m 

o 


j = 1,2, ... m - (3.27) 

i = n+4F +2m+l, . . . ,N+4N +3m 


th 

N = the number of transmission ratios from the j 
group , 

X. = the group characteristics. 


3.8 Solution Scheme: 

The object is to minimise f (X) subject to the N 

c 

constraints. The constraint list is far from being exhaustive. 
Even with these limited constraints, the solution becomes a 
formidable task because : 

a) The number of teeth on each gear must be an 
integer i.e. it could be 34 or 35 but never 
34.67. 

b) The availability of cutters restrict the usage 
of standard modules. This is a situation where 
some of the variables can only take predeter- 


mined values 
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These conditions make the formulation that of a 
Nonlinear Prctlem with Integer and Real Variables- Some 
solutions for pure linear integer programming problems have 
been obtained but the nonlinear part is still unexplored. 

^he only feasible remedy for this impasse seems in 
floating the variables. After reaching the minima by any of 
the existing techniques the variables are given their next 
higher integer values. Eut this method Is also not free of 
faults. The final variables must satisfy all the constraints. 
Depending on whether the module value is made higher or lower 
to achieve the next standard value, the number of teeth of 
corresponding gears should be lowered or raised respectively. 
This transformation of variables also requires that the com- 
patibility criterion is satisfied (i.e. the total number of 
teeth for any two pairs of mating gears between the two 
shafts must be the same). Sven now, the possible combination 
which yields the minimum centre distance may not be acceptable. 
The output speed at the spindle must not differ by more than 
two par cent from the required speod. Since the rounding-off 
of the teeth yields a change in the transmission ratio, there 
is a consequent change in the output speeds. But these may 
not be as difficult as they appear to be. Because of the 
finite number of design variables, the simplest method coulcL 
be to check every possible combination by fixing the range 
for each variable around the mean value found and always 
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checking all the constraints. That ’which leads to minimum 
centre distance amongst the acceptable combinations is the 
best speed diagram arrangement. 

3.9 Penalty Function: 

The problem is converted to one of the unconstrained 
minimisations by generating a function P (X) of the design 
variables such that 

N 
' c 

P (X) = f (X) + r t: \ (3.28) 

i=i %ao 

At any constraint surface i, (X) reaches zero thereby 

shooting up the values of total function P (X). The function 

30 

P (X) is given as high value z s 10 in case any of the cons- 
traints is violated. This helps in identifying the situation 
that one or more constraints are violated. In the algorithm, 
the solution then starts from the previous feasible design. 


/ 












CHAPTER 4 


THE SOLUTION ALGORITHM 


The technique employed in the solution of the 
foregoing constrained minimisation problem is discussed in 
this chapter. The constrained problem is first converted to 
an unconstrained minimisation problem and is then solved as 
a sequence of unconstrained problems. 

4.1 Conditions for Minimum s 


Necessary and sufficient conditions for a minimum 
of a function f (X) of n variables is that there exists X* 
for which 

f (X*) < f (X) (4.1) 

for all X. If the inequality holds in some neighbourhood of 

— ★ 

X then the function is said to hove a local or relative 
minima at that point. 


are 


The necessary conditions for the relative minimum 


3f 
3*i *. 


= o 


i = 1, 2, 


(4.2) 


The set of nonlinear simultaneous equations, thus obtained, 
are not easy to solve. Further, even if the above set could 
be solved there would be no guarantee that the obtained 
solution is not a maximum or saddle point rather than a minimum 
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One way to overcome this difficulty is to reduce the cons- 
trained minimisation problem to an unconstrained minimisation 
problem, 

4.2 Conversion to Unconstrained Minimisation : 

A constrained optimisation problem can be cast as 
a Mathematical Programming problem in the form : 

min f (X) 

subject to 

g. (X) ^ 0, 3 = 1, 2, ... K e 

where X is an n-dimensional vector of design 1 variables 
X. , i = 1, 2, ... n and g (X) are the given constraints 

3 

on the design. The function f (X) is called the objective 
function and its choice is governed by the nature of the 
problem. 

The constrained minimisation problem is converted 
to an unconstrained minimisation problem as follows : 

N 

P ( X, r ) = f (X) + r H — 1— (4.4) 

i=l g ± (X) 

where P (X,r) is the penalty function and r is an arbitrary 
penalty parameter which in the limit goes to zero. 

The minimisation proceeds over a strictly monoto- 
nically decreasing sequence of r - values from an initial 
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design vector X Q essentially inside the boundary of the 
constraints and for a positive value of r, say r^* The 
minimum of P (X, r ) must lie inside the constraint set 
because at the boundary some of the (X) tends to zero and 
consequently P (X) tends to infinity. The minimum of 

P (X, r ) depends on r , the starting value of r, and can 

1 

be written as X (r^). By reducing r^ in the next iteration, 
the summation term which penalises for the closeness to the 
constraint boundaries is reduced and hence in minimising 
P (X, r) more emphasis is placed on f (X). The method has 
special advantages of not following the zig-zag pattern of 
minimisation (Pig. 4.1) as is observed with other nonlinear 
p?~ogramming algorithms. This is because the minimisation 
starts from within and does not follow the boundary. 

( 12 ) 

Fiacco and McCormick have shown that if 

(1) the interior of the constraint set is non empty, 

(2) the functions f (X)and (X), i = 1, . . *N C 

are twice continuously differentiable, 

(3) the set of points in the constraint set for 
which f (X) v Q is bounded for every finite 
v , and 

(4) the function f (X) is bounded below for X 
in the constraint' set 

then, the optimal solution to the unconstrained 
problem approaches a local minimum of the 
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constrained problem as the value of r approa- 
ches zero. If, in addition, 

(5) f (X) and - g. (X), i 1, 2, N are 

i c 

convex functions, and 

(6) P (X , r ) is strictly convex in the interior 
of the constraint set for every r 

then the optimum solution to the unconstrained minimisation 
problem approaches the absolute minimum of the constrained 
problem as r approaches zero. 

4.3 Choice of Initial r : 

To start the algorithm the first decision is to be 
taken about the value of r and the factor C, the factor by 
which the values of r shall be decreased. The algorithm 
imposes no restriction on the values of r and C except that 
r y 0 and C 1. With large r - values the number of compu- 
tations required to reach the constrained minimum increases 
tremendously while with small value of r, the function is 
highly distorted and requires initial solution very close to 

, _ . . Cn ) 

actual minima. 

Two methods to choose the initial value of r have 
been suggested in the literature. , 
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T 


( 1 ) 


where 


f (x 0 ) . p (x 0 ) 


(4.5) 


N. 


« ) - H 


p a 0 ) 


1.0 


5=1 


g i 


and represents the gradient of the function. 


Case I : r 0 
1 


If r 0, minimisation of f (X ) alone, without 
1 >• o 

considering the penalty term, is carried out. At every new 
point the r-value is checked by equation 4.5. If the value 
is positive, unconst rained minimisation with obtained velue 
of r is carried out. 


Case II i r = 0 
1 


This means that the unconstrained minimum has been 
_ 

reached and X = X 
o 


Case III : r^/’O 


If r > 0, this is taken as the starting value for 
usual minimisation. 


( 2 ) 


That value of is selected which makes the objec- 
tive function and penalty term equal that is, ** 


r 

1 


f « 0 ) 



1.0 

gl (30 


(4.6) 
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This has two weaknesses. Firstly, the starting 


point may be too close to some of the boundaries thereby 
making such selection of r useless. Secondly, there may not 
even be a point in the feasible domain which is not near a 
constraint boundary. If some of the constraints are near the 
boundarjr, i.e. penalty is large the value of r obtained from 
equation (4.6) is very small. This difficulty is eliminated 
if in the summation for penalty those constraints which are 
quite close to the boundary are substituted by some minimum 
value of constraint. Modified expression for r^ is 



(4.7) 


Since the constraints are normalised here, the value of g . 

5 to min 

can be arbitrarily selected as 0.2. 


4.4 Extrapolation: 


The vectors X (r^), X (rg) ..., X ( r^_) obtained by 
minimising P (X, r) over decreasing sequence of r - values 
lie on a trajectory X (r), where X (0) is the desired solution. 
From this trajectory, it is possible to obtain, by extrapola- 
tion, estimates of final solution X (0) and the next minimum 
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It has been shown that the trajectory approximates 

i ( 8 ) 

The polynomial of r . This is because as r tends to zero, 

JL 

the function behaves as a linear function of r 2 . The appro- 
ximating function is defined as 


jl (k-l)/2 

X (r) = a Q + a^ r* + a- 2 r + . r 


(4,8) 


where 


a^’s are the undetermined coefficients. 


For small v. lues of r 


X (r) X (0) + a 
Solving these for X (0) 



(4.9) 


X (0) 


X (r/c) - X (r) 



(4.10) 


An estimate of minimum point for next r - value is obtained 
by assuming that 

X (r/C 2 ) ' X (0) + a (4.ll) 

Solving equations (4.10) and (4.11) for X (r/C 2 ) 


X (r/C 2 ) = X (r/C) + X (r/C) -X (r) (4.12) 

This value of X can be used for the next minimisation, thereby 
substantially reducing the effort for minimisation of P (X), 
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However, these extrapolated values have to he first checked 
so that they do not violate the constraints. Estimates of 
X (0) he coxae reasonably good after three or four trials, and, 
towards the end offer more accurate estimation of values of 
X (0). 


4.5 method of Unconstrained Minimisation : 

(Variable Metric method of Fletcher and Powell) 

The method is based on the properties of a quadra- 
tic function and is designed so that when applied to a 
quadratic, it minimises the function (of n variables) in 
n iterations. 


Central to the method is a symmetric positive 
definite matrix M which is updated at each iteration. It 
collects in itself the information about the curvature of 
the function P (X) and supplies the current direction of 
move d^ by multiplying it with current gradient vector. 

f A typical iteration in the minimisation of f (X) 
proceeds as follows : 


(l) given the starting point iL and the gradient of 
P (X) at X^, d^ (X), the direction of (I + l)th 
minimisation of P (X) is given by 



( 4 • 1c ) 
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(2) find <X[* so that f ( Zj_ + <2[* d^ ) is minimum 


along the line d^. 


(3) set .X 


i+1 


X i + °1* d i 


(4) calculate the new gradient vector 

Vi =VPd W 


and set 

Y = G. 
i i+1 


Lr . 


(5) • calculate the new j^Hj| matrix, i*i] by 

r 1 r 1 d. \ tel Y Y^ tel 

[ H i [ H i] * .LiL , 1 lLiJ 


d. Y. 


w 


(4.14) 


(6) begin the next iteration from (l). 


It can be shown in the case of a quadratic function 
that matrix £ converges to the inverse of the matrix of 
second partial derivatives of the quadratic after n iterations. 
In case of a general function, M tends to the inverse of 
the matrix of second partials of the function evaluated at 
the minimum. This is because as the minimum is approached 
the second order term in Taylor series predominates. 


If during any iteration the minimum is not found 
along d^ then is not modified in the next iteration. The 

new direction is calculated from equation (4.13) using 
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and at any step jjsQ is no longer found to be 

positive definite then it is reset to unity. The positive 
(13 ) 

definiteness at point X. is tested by evaluating 
T 1 

Xi X r If this is positive then is positive definite. 


4,5 The Gradient: 


At every iteration the gradient has' to be calcula- 
ted in the minimisation algorithm. It can best be determined 
by -writing down the exact derivatives of P (X) with respect 
to the design variables. Eut , as in the present problem, 
both the number of design variables n and the number of cons- 
traints depend on the problem, one cannot arrive at a general 
analytical expression for the gradients. Moreover even in a 
moderately complex engineering problem, the analytical gra- 
dient calculation is very cumbersome. Gradients are there- 
fore approximated by the difference technique. 

The gradient of the function P (X) with respect to 
the design variables is given by 


G (X) 


VP (x) = ( 




a p ) T 
>ax n 


(4.15) 


I. j, J. ; - r r\m 

CEN1KM. ! iffRARV 

T 54874* 


No. 


It may also be written as 
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p (X) = c - r 2 , ... 

i=i g 

i 


ft. 


2-n 


y* JL 's 

Z_J 2 • 0X J 

■i =t 6 . £*• 


(4.16) 


i=l 


Expanding the function F (X) of one variable about 


X by Taylor Series 
o 


d£ 


'1' * '""o' ’ ' dX 


F (X J = F (X J + ( ) 

. 2 . 


X Q - V V + 


i ( 


d F 
dX 2 


) 


x d C X 1 ~ X 0 } + *** 


(4.17) 


neglecting the higher order terms 


F (X n ) - F (X n ) 

dF = 1 0 (4.18) 

®lx 0 - x 0 > 

Close to the minimum, the function P (X) has very 
large curvature indicating significant values for the second 
and third derivatives which have been neglected in the forward 
difference scheme. Even if g (X) and f (X) are nonlinear, 
still they are quite well behaved as compared to P (X) at any 
point inside, the domain. Thus, the gradient from the equation 
4.16 should be used. 


4.7 Linear Minimisation 

In step 2 of the Fletcher and Powell method it is 
essential to find S £J* for which P ( X^ + <5* d^ ) is 


minimum 
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If X^ and d^ at particular i are taken to be cons- 
tant; then P ( + (F d^) is a function of CT alone. 

p ( x ± + cr d i ) = f(<r) (4.i9) 

Approximating the function ^ ( 0~) by a cubic as shown in 
Fig* 4.2 the minimum of ^ ( CT ) can be approximated by 0~* in 
the following way. 

1. The initial step size is computed as 

f " 2 F P a i ) ’ 

CT£ = Min j 1.0, 

L y ( (T) . d. 

where $ is the lower estimate of the function 
min 

f (X) . It is usually assumed that in each minimisa- 
tion the function value will reduce by 5 to 10 percent. 

2 . Evaluate 

T 

ty ' - d^ VP ( X^ + d^) at the points, 

°^i = OjCr^ > 2<3 i» 4c5 i> ..., a , b. 

T 

where b is the first of these values at which ^ 
becomes non-negative or ^ has not decreased and a 
is the point preceding b. 

3. A cubic polynomial^ ^is now fitted to the four 
values of ^(a), W (a), ^(b), and ty* (b). Its 
minimum is taken as 0~ . 


] 


(4.20) 
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. (b - a) 

(4.21) 

whe re 

z = 3 r + y/ u) + f' (b) 

r — x 

¥ = Z 2 - f' (a) . (b) 

If (Ch) is less than t^(a) and (b), (7^ is accepted as 

<T* . Otherwise, depending on whether ^ ( O' ) is positive or 

c 

negative, the interpolation is repeated over the sub interval 
(a,cr e ) or ( G~ e1 b) respectively. 

At 0^=0*, i{j' (CT) =0 which means that the 
direction of the new gradient is orthogonal to the previous 
direction, ^ven if the check 


-J^±jU L=Z 

*Kb) - £Ta) + 2 ¥ 



4.8 Slope in d Directions 

In one dimensional minimisation of P (X), slope 
in the direction d is required. This can be 'Obtained in two 


ways: 
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!• By taking dot product of the vector d and G at 
the point 

hi . 1 — ?■ — 

Y = d . G (4,22) 

2. By taking finite difference in d direction. 

In the first method, to calculate the gradient, 
it is essential to evaluate the function value n times, 

Where as, in the second method only one computation of the 
function is necessary for the forward difference scheme. 
Therefore, the second method has "been preferred here. 


The forward difference scheme in the d direction is 


as follows 


H = x ± - e.d ± 


where ^ is a small step in d direction. 


. r — 

The slope y ( (T) in d direction is then 


N 


{£/' _ d ¥ (O') _ d f (X) ^ V 

r - - “a? " r Z_j 


Hi 

d£~ 


i=l 


g 


( 4 . 23 ) 


where 


.df (x) = fjx!l_-_xxxl 

d£ € 

d g i g ± <X*) - gj_ CX) 

d€ " £ 
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4.9 Convergence Criteria: 

The following are the convergence criteria used 
in the algorithm at various stages. 


1 . 


For one dimensional minimisation 
I T 1 


a) 


a . o_ _ 

idljCh 



(4.24) 


where 



is the cosine of angle between d and G. 


h) the percent change in the function is less than 
specified tolerance 

I f ub i < * ( 4 - 25 > 

I f CD 1 ^ s 2 

2. For minimisation for a particular r 

a) when the approximation for the percentage change 
in P (X) in the next iteration has become less 
than a certain pre-as signed value 


g C h3 g 

P (X) 



(4.26) 


b) another criterion to stop the iterations for a 
particular r is when n/2 consecutive one dimen- 
sional minimisations have yielded CT (step in 
d direction) less than one percent of the sma- 
llest value of the design variables . 
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3. 


For the constrained minimum 

a) when the percentage change in the objective 
function between minima corresponding to two 
consecutive r values is less than the given 
tolerance, that is 


mm 


(X 


r-l) - 
^min 


f min 

TFT 


( x 




2 


(4.27) 


b) when the minimum percentage change in any design 
variable corresponding to two consecutive minimum 
r values is less than the given tolerance 


*i ( r . C ) - X ± ( r ) 

» i 

c) when the percentage change in the extrapolated 
function value f (0) is less than the given 
tolerance . 

f ( r ) 


< ^r 


(4.29) 


< £ (4.28) 

\ rn 






MINIMIZATION 


FIG' 4 2 


CHAPTER V 


RESULTS AND DISCUSSIONS 

Nine illustrative examples have been considered 
to show the applicability of the constrained minimisation 
algorithm to solve gear box problems. First example is a 
main drive of an actual lathe. Results obtained therefrom 
are tabulated. Other examples have been studied to select 
the best layout diagram for an 18 - speed gear box. The 
criterion for "best" design in all the discussed examples 
is the minimum centre distance. Constants which are common 
to all examples are tabulated in Table 6.1. 

5.1 Example r 1: 

This is the main drive of a 10-speed gear box 
as shown In Fig. 1.3. It has 5 groups and also bach gear 
arrangements with composite gears. 1st, 3rd and 5th group 
can have one independent driver gear therefore first six 
design variables denote number of teeth and their modules 
alternately. 7th, 8th and 9th variables indicate indepen- 
dent transmission ratios and account for the fixation of 
speed diagram. Input i nformat ion for this problem is given 
in Table 6.2. Table 6.3 compares the existing and opti- 
mum values of gear teeth, their modules, minimum transmi- 
ssion ratio from each group and individual centre distance 
of each group. 



the centre distances are shown in Table 5.6. Initial and 
final values of lowest transmission ratios are also indi- 
cated. The usual time for one complete run of each exampl 
is between 12 to 14 minutes. Since the component designs 
in these cases will not provide any useful information 
their computations have been avoided to save computer time 
The examples have the same input information as example 
2 except for the quantities indicated in Table 5.7. Ini- 
tial and final speed diagrams for all the problems are 
indicated in Figures 5.5, 5.6, 5.7, 5.8 and 5.9 respec- 
tively. 

5.4 Example - 8: 

This is the same gear box as in example 7 but 
with different starting design vector. The initial and 
final design variables are indicated in table 5.7. 

5.5 Example - 9: 

Same example 7 is solved for 100 rated horse 
power capacity. The results obtained necessarily have 
greater centre distance as compared to example 7 but the 
final speed diagram does not change considerably from the 
optimised results of speed diagram of example 7. Optimum 
point has 7 bounded constraints. 



5.6 Discussions : 


The results from 1st problem show about 16$. re- 
duction over centre distance from the existing design. The 
rounded number of teeth have adjusted themselves in such 
a fashion that out of two mating gears the number of their 
teeth have no simple common factor which is desirable 
from the woar considerations. This added advantage shall 
be obtained almost always because the transmission ratios 
from different groups do not arrange themselves in simple 
fractions. The final speeds ere within the tolerance -of 
the preferred speeds. Centre distance in the first group 
in example 1 remains almost the same because the speed of 
1st shaft does not change. The effect of change in speed 
diagram is observed in the later groups. Study of initial 
and final speed diagrams of different examples clearly show 
that the intermediate shaft try to attain as large speed as 
possible. This may be the cause that example 7 and 8 reach 
the same speed diagram. Results from example 9 indicate 
that there exists only one speed diagram for a particular 
layout which shall always be reached for any horse rower 
capacity and same input and output speeds. Results from 
example 3, 4, 5, 6 and 7 indicate that the layout diagram 
used in example 5 for 18 - speed gear box is desirable 
from considerations of minimum centre distance. 
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The percentage improvement on the objective func- 
tion in the first example may be attributed to one or more 
of the following points. 

1. The gear box was designed to suite the re- 
quirements of already existing castings. 

2. The method of gear design may be different 
from the one used here. 

3. The design was not optimum oven for the 
chosen speed diagram. 

4. No attempt was made to exploit the independent 
design parameter i.e. transmission ratio to 
aid in the reduced size design. 

The optimum search is possible here because the 
algorithm seeks the minima (at least relative) in the fea- 
sible design space. 

The author does not claim that a global optimum 
has been reached for example 1. This is because firstly 
not many starting vectors were chosen due to shortage of 
computer time and secondly even for one starting point the 
iterations were stopped when no significant changes occured 
in objective function. Actually this was the trade off 
between computer time and accuracy desired. 
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TABL'D - 6.1 


VALUES OF C OF STATUS 

COMmOTI TO ALL EXAMPLES 

Identifier 

Magnitude 

SIG W 

20 Kg/mm 2 

SIG C 

12 Kg/mm 2 

SIG T 

8 Kg/iorn^ 

FSC 

0.8 

FNUT 

1.2 

FS 

1.25 

E 

2100.0 Kg/mm 2 

SC 

1.125 

ALH 

10.000 hrs. 

HOUR 

4.00 

ANG ( ) 

180° 

GDCH 

1.3 

GDP 

1.7 

RHO 

0.000008 Kg/mm 2 

PLF 

175 Kg. 

CUTF 

200 Kg. 

DSP 

32.0 mm. 

DSJ 

35.00 mm. 


AAMOD ( ) 


SPSl 
BPS2 
EPS3 
EPS4 
IT IB. 


1.0. 1.125.1.375.1.5.1.75.2.0.2.25, 

2.5.2.75.3.0. 3.25.3.5.3.75.4.0, 

4. 5.5. 0. 5. 5. 6. 0.6. 5. 7. 0.7. 5. 8.0, 

9.0. 10.0, 


0.005 

0.001 

0.100 

0.005 

100 



TABL‘D 6.2 


VALUES OF THE CONSTANTS USED IN EXAMPLE I 


Identifier 




Values 



HP 




12. 

0 



NSH 




6 




ngmax 




5 




IEGSH 




4 




i lap 




4 




IP ( ) 


(3, 2, 

If 1 

i 1) 

T 


ICAR ( ) 


(1, 2, 

1» 1 

f 1) 

T 


ICOi'iP ( ) 


r-** 

0 

0 

0 

0 

0 




0 

2 

0 

2 

0 




0 

0 

0 

0 

0 




0 

0 

0 

0 

0 




0 

0 

0 

0 

0 




_0 

0 

0 

0 

0 


DIS ( , ) 

325 

160 

135 

1 




230 

140 

307 

140 

120 


280 

27 

140 

1 


1 


. is 

0 

-33 

1 

1 


1 


150 

75 

1 

1 


1 


450 

1 

1 

1 


1 
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TAbLE 6.3 


THE TABLE COMPARES THE 

i 

& 

M 

K 

O 

DESIGN TO 

THE EXISTING ONE 



Existing 

Optimum 



De sign 

Design 

1. Number of teeth/module 

S.No. 

1,1 

1,2 

25/3.0 

31. 11/2.73 


30/3.0 

36.92/2.73 


1,3 

2,1 

35/3.0 

42.55/2.73 


15/3.0 

18.33/2.73 


2,2 

30/3.0 

30.70/2.73 


2,3 

35/3.0 

36.51/2.73 


2,4 

30/3.0 

30.71/2.73 


2,5 

25/3.0 

25.07/2.73 


o ’ "1 
° , J- 

25/3.0 

35.31/2.21 

, 

3.2 

3.3 

53/3.0 

30.58/2.73 


38/3.0 

18.21/2.73 


4,1 

25/3.0 

25.67/2.21 


4,2 

50/3 . 0 

55.35/2.21 


5.1 

5.2 

47/3.5 

32.34/2.25 


50/3.0 

64.99/2.21 


6,1 

47/3.5 

84.75/2.25 

2. Lowest transmiss ion 

1 

0.717 

0.852 

ratio from a group 

2 

0.283 

0.59S 

o 

O 

0.500 

0.638 


4 

0.500 

0.395 


5 

1.000 

0.386 

3. Centre distance of 

1 

90.0 

92.3 

groups 

2 

99.0 

66.5 


O 

W-' 

112.5 

96.0 


4 

112.5 

96.0 


5 

164.0 

151.5 

L. Total centre distance 


578.5 

482.3 
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TABLE 6.4 

ROUNDED OFF GEARS AlTD MODULES AND DESIGN OF COMPONENTS 

A 


Gear 3 
leeth 3 

__ L 

Modules 3 

5 

(5 

Corrections ^Crushing 5 Z * 3 Limits 

0 Area 5 m 0 Nominal IT Lower ~ 
3 3 3 Value 3 

3 Upper 

31 

2.75 


0.00 


15.6 

4 

25.05 

24.99 

25.01 

37 

2.75 


0.00 


15.0 

5 

37.80 

37.74 

37.75 

43 

2.75 


-0.23 


15.0 

5 

37.60 

37.54 

37.56 

18 

2.75 


+0.33 


16.6 

3 

21.55 

21.49 

21.51 

31 

2.75 


0.00 


16.6 

4 

25.05 

24.99 

25.01 

37 

2.75 


0.00 


15.0 

5 

37.80 

37.74 

37.76 

31 

2.75 


0.00 


16.6 

4 

25.05 

24.99 

25.01 

26 

2.75 


-0.15 


15.0 

3 

20.95 

20.89 

20.91 

36 

2.25 


0.00 


24.0 

5 

30.09 

30.03 

30.06 

31 

2.75 


+0.18 


16.0 

4 

29.80 

29.74 

29.76 

19 

2.75 


0.00 


16.0 

vS 

20.95 

20.89 

20.91 

26 

2.25 


0.00 


2S.0 

O 

O 

17 .32 

17,26 

17.28 

55 

2.25 


0.00 


29.0 

7 

44.80 

44.74 

44.76 

33 

2.25 


0.00 


54.0 

4 

24.20 

24.14 

24.15 

65 

2.25 


0.00 


54.0 

8 

51.60 

51.54 

51. 56 

83 

2.25 


0.00 


57.0 

9 

59.00 

58.94 

58.95 

; Number of 

teeth over which 

measur 

ement 

is to he 

taken . 



B 

Group No. 3 Centre distance 3 Lowest Transmission 
Lxatio 


93.50 

0.837 

68.75 

0.580 

102.50 

0.655 

102.50 

0.400 

130.50 

0.398 


5 
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TABLE 6.4 

ROUNDED OFF GEARS AND MODULES AND DESIGN OF COMPONENTS 

A 


f GearT Modules $ Corrections §C rushing 5 Z * § Lim its 

5reeth8 5 8 Area 0 111 5 Nominal 5" Lower (T Upper 

XX X Y Y Y V Y X X 


1 

y 0 

1 



5 Value 

JL 

JL 

31 

2.75 

0.00 

15.6 



4 

25.05 

24.99 

25.01 

37 

2.75 

0.00 

15.0 

5 

37.80 

37.74 

37.75 

43 

2.75 

-0.23 

15.0 

5 

37.60 

37.54 

37.56 

18 

2.75 

+0.33 

16.6 

3 

21.55 

21.49 

21.51 

31 

2.75 

0.00 

16.6 

4 

25.05 

24.99 

25.01 

37 

2.75 

0.00 

15.0 

5 

37. SO 

37.74 

37.76 

31 

2.75 

0.00 

16.6 

4 

25.05 

24.99 

25.01 

26 

2.75 

-0.15 

15.0 

3 

20.95 

20.89 

20.91 

36 

2.25 

0.00 

24.0 

5 

30.09 

30.03 

30.06 

31 

2.75 

+0.18 

16.0 

4 

29.80 

29.74 

29.76 

19 

2.75 

0.00 

16.0 

r ^ 

C 

20.95 

20.89 

20.91 

26 

2.25 

0.00 

29.0 

O 

17 .32 

17.26 

17.28 

55 

2.25 

0.00 

29.0 

7 

44.80 

44.74 

44.76 

33 

2.25 

0.00 

54.0 

4 

24.20 

24.14 

24.15 

65 

2.25 

0.00 

54.0 

8 

51.60 

51.54 

51 . 56 

83 

2.25 

0.00 

57.0 

9 

59.00 

58.94 

58.95 

: Number of teeth 

over which 

measur 

ament 

is to he 

taken . 



B 

Group No. 5 Centre distance!) Lowest Transmission 

... L_ 5. ratio 


1 

93.50 

0 . 837 

2 

68.75 

0.580 

3 

102.50 

0.655 

4 

102.50 

0.400 

5 

130.50 

0.398 



TABLE 6.4 (CONTINUED) 
C 
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S .No .~F . UHAFT_. EBAR lflGS ~ g ~ C LUTC EB's 



5 Dia- 

1 Deflec- 

1 Static 

5 Dynamic g 

Type 

g T hernia 


fi meter y tion 

X X 

LCapa_citx_ , LJL^acitv , j 


g Utilis 


_JL 

JL 

-LL.eft ±_ 

JligfcjL LLeXtJ Jilgilti 




1 

32 

0.039 

504 

455 2185 

1207 

SK-lOd 

45.76 

2 

32 

0.025 

863 

480 2040 

1200 

EK-lOd 

59.82 

3 

35 

0.017 

785 

783 2224 

1979 

BK-lOd 

93.31 

4 

36 

0.031 

708 

1424 2283 

3593 

SK-20dD 

28.61 

5 

38 

0.007 

917 

1500 2170 

4500 

EIC-20d 

98.83 

6 

50 

0.032 

165 

2000 504 

6100 




D 




Preferred 

Speed g Actual, Speed 

JL_ JL E.r ror. 

1 

45 

46.0 

+2.20 

2 

64 

65.7 

+2 . 63 

3 

90 

92. 0 

+2.21 

4 

127 

130.0 

+2.30 

5 

180 

177.0 

-1.70 

6 

255 

253.0 

-0.80 

7 

355 

350.0 

-1.40 

8 

500 

497.0 

-0.60 

9 

710 

710.0 

% 

0.00 

10 

1000 

983.0 

-1.70 



TABLE 6.5 

INPUT DATA FOR BXAI-iPLE 2 
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Iden tifier 


Value 


RPMI 

RPML 

RPtoH 

HP 

NSH 

N GMAX 

IBGSH 

ILAP 

IP ( ) 

ICAR ( ) 

ICOmP ( , 


0 

0 

0 

0 

0 


1400.0 

35.5 


1600.0 

12.0 

5 

6 
0 
0 

(3, 3, 1, 
(1, 3, 1, 

0 0 0 
0 0 0 
0 0 0 
2 2 0 
0 0 0 



0 

0 

0 

0 

0 


0 

0 

0 

0 

0 


TABLE 6.6 


COMPARISON OF EXAmPLE NO. 3, 4, 5, 6 c 7 


Sjj 

No. 

Ex. 3 

Ex. 4 

Ex. 5 

Ex. 6 

Ex 7, 

..ICAR ( ) 

1 

6 

2 

1 

1 



2 

2 

6 

3 

6 

1 


3 

1 

1 

1 

1 ■ 

1 


4 

1 

1 

9 

3 

*->1 1 

O 

. Initial 

1 

0.30 

0.80 

0.80 

0.30 

0.30 

U ( ,1) 

2 

0.80 

0.30 

0.50 

0.30 

0.80; 


3 

1.00 

0.50 

0.50' 

1.00 

1.00 


4 

0.45 

0.8S 

0.81 

1.17 

0.35; 

s. Final 

1 

0.37 

0,76 

0.89 

0.90 

0.36! 

U ( ,1) 

2 

0.68 

0 .36 

0.71 

0.36 

0 • 8«3; 


o. 

1.10 

1.05 

0.55 

0.97 

0.97 


4 

0.39 

0.37 

0.47 

0.34 

0.36' 

1. Centre 
Distance 


362 

314 

303 

320 

351 ! 

.:nti oium.ui .mm. t i,. ir , *,•«» * 

w , , « 

^ . — — 

4 

- • • 


* - - *’ * ~ 



TABLE 6.7 

COMPARISON OF EXAMPLES 3 AND 8 
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EXAMPLE 3 EXAMPLE 8 



Initial 

Final 

Initial 

Final 

Design 

Variable 

0.161 

0.069 

0.563 

0.080 

(Normalised) 

0.889 

0.122 

0.421 

0.420 


0.161 

0.122 

0.333 

0.127 


0.889 . 

0.204 

0.421 

0.329 


0.161 

0.098 

0.213 

0.110 


0.889 

0.220 

0.617 

0.230 


0.028 

0.064 

0.302 

0.081 


0.314 

0.334 

0.444 

0.33S 


0.420 

0.413 

0.342 

0.409 






OPTIMAL Speed J)ia^ram For. Example No: 1 
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Fi<5. 5-2 






ooti 
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INPiJT Speed Diagram for Example Mo-.2 

r i g . , „ 




















CHAPTER VI 

CONCLUSIONS AND RECOMMENDATIONS 

The application of the proposed method to solve a 
gear box problem shows that the complicated machine design 
problem can be tackled by computer. That the algorithm is 
of any practical use or not depends on the problem and the 
situation in which it is encountered. In general} the ques- 
tion as to whether a mathematical programming approach is 
suitable or otherwise for a practical design problem depends 
on the following considerations : 

1. The need for obtaining an optimum design. 

2. The effort and cost required to write a pro- 
gramme identifying the variables and cons- 
traints in the problem compared to man hour 
requirements in a conventional design office 
for the same problem. 

3. The capability of the programme to tackle a 
wide range of similar problem. 

The programme that has been developed is very 
general and can handle any speed gear box of conventional 
lay out diagrams. Therefore it spreads out the cost of de- 
veloping the programme and the computer time required to 
achieve a satisfactory run. Moreover the design office re- 
quirements to completely design a gear box is tremendous. 
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The author does not claim to give sufficient information for 
final manufacturing of the gear box, hut then, the results 
of the programme defenitely allow the designer to draw the 
final drawing and obtain other important dimensions from 
the assembly drawing itself. Thus the routine part of the 
analysis Is not needed and using computer the designer ob- 
tains a optimum design as a bargain. Finally, it can be 
stated that the time and cost of development of the general 
programme, though considerable is quite insignificant to 
that for preliminary design of a single gear box. Many po- 
ssible layout diagrams for the same speed can be studied 
before accepting one of them. 

6.1 Recommendations : 

The author recommends the following as a basis for 
further work in the field: 

1. One very severe short coming of the algorithm 
is that it is not self starting. One has to 
feed starting feasible solution. However, 
remedying this is not as difficult as it sounds. 
The initial number of teeth and module can be 
given as very high values thereby making no 
room for inadequate power transmission. Only 
difficulty remains in establishing initial 
value of independent transmission ratios. 
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Generally putting the values between 0.5 to 
1.0 in any order has resulted in a feasible 
design but it may not be true universally. 
Dependent transmission ratios have to be cal- 
culated to 'check the constraints on them. 

Had it been a self starting algorithm, this 
trouble would have been avoided, 

2. The author strongly feels the necessity of 
main machine interaction at some of the' 
places in the gear box design. One example 
is fixing of distances of the gears from the 
left hand bearing. This has been done here 
on the approximate judgement of the designer. 
However, , the human factor can be dispensed 
with if variables along the axis of the shaft 
are introduced. These will fix the Z coordi- 
nates of all the gears along the shaft with 
clearance as instructed,' 

3. The programme shall be more useful if the 
angular disposition of shaft is also intro- 
duced as a variable in space optimisation. 

But this will lead to other assembly restric- 
tions. One of them is the location of the 
guide ways for movement of gear change levers. 
Exact assembly restrictions be incorporated 
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in consultation with Machine Tool manufactu- 
ring concerns. 

4. In gradient calculations, the derivatives 

of constraints which do not vary with respect 
to certain design parameters have also to he 
computed for uniformity, though it is obvious 
that the derivatives would be zero. This 
accounts for much computer time. An algorithm 
may be developed in a minimisation programme 
which would identify those constraints which 
do not vary with respect to a variable about 
which the derivative is being calculated and 
skips the derivative computation. 
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APPMDIX A 


In this Appendix it is mathematically proved that 
the three composite gears as a shaft are not possible. The 
maximum and minimum output speed with two composite gears are 
found. This analysis is to stress the use of fewer pairs of 
gears for the same transmission. 

A.l Impossibility of Three Composite Gears on a Shaft: 

To arrive at 18 spindle speeds the minimum numbers 
of gears is shown to be 11 in the Fig. (A.l). The problem 
which has been discussed earlier used 17 gears for the same 
purpose. But not all the rest are redundant or disponsible. 

It has been shown below for the arrangement in Fig. (A.l) that 
the combination will lead to active number of teeth for few 
of the gears if we fix one of the gear teeth say a. In; fig. 
the different gear ratios are indicated. S is the minimum 
desired speed and the input speed is taken to be 1 for conve- 
nience. The other speeds on the shaft indicates the fractions 
of the input speed. The letters a, b, c etc. are used for 
actual number of teeth on those gears. 

Layout diagram for the first two groups in the 
Fig. A.l is shown in Fig. A.la. Nine independent equations 
derived therefrom are indicated below: 
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s 

= 

£LslS. 

d.h 

(A.l) 

X 1 

s 

= 

b 

h 

(A.2) 

2x-, 

s - 1 

as 

C ,e 

f .h 

(A. 3) 

X 2 

S 6 

55 

a 

g 

( A«4) 

2x 0 

8 ^ 

= 

af 

di 

(A. 5) 

a + d 

= 

b + e 

(A. 6) 

a + d 

=s 

e + f 

(A.7) 

d + g 


e + h 

(A. 8) 

d + g 

- 

f + i 

(A.9) 


These constitute nine equations in nine gear teeth 
unknowns. Therefore, none of the gears can he arbitrari-ly 
chosen* By solving these equations, one gets a relation bet- 
ween a and d, which is not admissible for any combination of 
S and 0. 

A. 2 Analysis With Two Composite Gears: 

Instead of tackling the whole problem we shall 
see if even six spindle speeds are possible from the arrange- 
ment shown in Fig. A.2 which uses 8 gears. If it is, then 
it can be used as a subassembly for an 18 speed gear box. In 
that case the shaft number 1 in Fig. A.2 will have three 
different speeds by a cluster gear set from the motor shaft, 
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Equations (A--.3 ) and (A-14 ) are redundant as they can he 
derived from other equations i.e. L.H.S. and B.H.S. match 
separately by the following relationship: 

(A-14) = U . 1S) 


Now we have seven independent equations in 8 un- 
knowns therefore all gear teeth are solvable in terms of "a” 
and the minimum speed desired S for a particular layout. 

The solution is lengthy and is carried out in the following 
steps. Eliminating c and b from (A-12) and (A-13) respectively 
and substituting in other equations, then eliminating f, e, 
g, h alternately and substituting in the remaining equations 
we get following relationship between a and d. 


a 


( 


x 2 
•? - 


2x-, 

0*1 


2x,+X 0 Xp x 9 3 

(S0 1 2 - S0 - 0 + 1) . 0 


(A-19) 


while the other teeth number are written below in terms of 
a and the known teeth number. 

, _ (a + d) ...a (a-20) 

h “ Xo x n 

S . 0 . (a0 1 + d) 


(a, + d ) »._a, — (A-21) 

x 2 x 2. 

S . 0 . (a + d0 ) 


X 1 + x 2 

- S , d . g . 0 / a 


e 


(A-22) 
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Equations (A-. '3 ) and (a- 14 ) are redundant as they can be 
derived from other equations i.e. L.H.S. and R.H.S. match 
separately by the following relationship: 

(A-14) = J^nl3 , ) . ^ l A , -lU ana ( A -15) 


Wow we have seven independent equations in 8 un- 
knowns therefore all gear teeth are solvable in terms of ,r a" 
and the minimum speed desired S for a particular layout. 

The solution is lengthy and is carried out in the following 
steps. Eliminating c and b from (A-12) and (A-13) respectively 
and substituting in other equations, then eliminating f, e, 
g, h alternately and substituting in the remaining equations 
we get following relationship between a and d. 


a 


( 


x 2 
i) - 


2x n 

0 ±1 


(S0 


2x^+x 2 


Xo Xq 

- S0 - 0 d 


(A-19) 


+ 1 ) . 0 


while the other teeth number are written below in terms of 
a and the known teeth number. 


h 


£a + d) .a 

Xo X-i 

S . 0 . (a0 1 + d) 


(A-20) 


= .ffl, + d ) (A-21) 

x 2 x 2 x 

S . 0 . (a + d0 ) 

X 1 + x 2 

= S . d . g . 0 / a 


e 


(A-22) 



f 


8S 


- s . a . h . /Va 

b = S.. g . 0* s 

2x-> 

c = S . h , 0 X 


(A-23) 

(A-24) 

(A-25) 


is evident that the arrangement is meaningful 
only with such S and 0 values which mate the R.fl.s. of U-19) 
positive . 


Particular Cases: 

(!) According to the arrangement shown in Fig, (A.2a) 
where x-^ = l, Xg = 3 we get from (A-19) 

a — 3 0 

4 3 2 • ’ (A-26) 

(S0 + S0 - 0 - 0 - l) 

It follows that 
4 3 2 

S0 +S0 -0 - 0 - 1 > 0 (A- 27) 

for 0 — 14, and for inequality (A— 27) to be true 
we must have 

5 3> 0.647 

to avoid excessive big gear in "d 11 the inequality 
should be sufficiently above zero. The m inimum 
of the maximum spindle speeds attainable theore- 
tically will be 0.647 x 1.4 5 = 3.72 times the 
input speed of shaft 1. 
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(ii) For the arrangement shown in Fig. (A. 2b) where 
X 1 = X 2 “ 1 we § e ' t the inequality 

2 32 S\0 

(0 + 1 ) . ( S0 +S0 + 1 ) X 0 > 

^ 0 >0 

{A. 28 ) 

which is satisfied for all values of S and 0. 
This is a very encouraging result. This shows 
that even if the recommended layout diagram is 
not followed the number of moving parts in the 
arrangement can be reduced by using this as a 
part drive for the gear box. 












APPENDIX B 


DESIGN OF COMPONENTS 

Manual design calculations for the components of 
the gear box is highly time consuming* Moreover, it becomes 
almost impossible to take into account all the conditions 
of loading and select the critical one. The situation be- 
comes even more complex with the increase in the number of 
spindle speeds. The design of some of the components sub- 
ject to the worst loading case is discussed in this Appendix. 

B.l Shaft Design : 

From assembly considerations, the shafts are usua- 
lly designed in the form of stepped cylindrical bars (Fig.Bl), 

In the design of shafts the following considera- 
tions are Important: 

1. The bearings supporting the shaft are considered 
as simple supports offering no moment. 

2. The driving and the driven forces are generally 
non-parallel and non-coplaher. Only in case of 
composite gears the driving and driven forces are 
coplanar. 

3. The number of loading conditions on a shaft is a 
product of the number of driver and driven gears 
on that shaft. In case of first shaft and the 
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spindle, it is equal to the number of driver and 
driven gears respectively. 

4. Since the least speed produces maximum torque and 
load on the shaft, minimum speed for particular 
loading condition of the shaft is considered for 
bending moment calculations. 

5. The resultant force on a gear acts along the pre- 
ssure line. 

6. Distances of gears on a shaft should be given in 
order to calculate the bending moment at different 
sections of the shaft. 

7. The relative disposition of shaft should be given 
as input information. In the present work this 
is given in terms of included angle in clockwise 
direction. Fig. B2 shows the included angles. 

8. A suitable coordinate system should be chosen to 
facilitate the calculation of the bending moment 
of shaft in two dimensions. A shaft can be loaded 
by only two forces at a time, namely the driver 
and the driven force. The co-ordinate system 
chosen here is an oblique system with X and 1 
axis coinciding with the direction of the two 


forces . 
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There may be three distinct types of loadings for 
a gear box shaft. 

a. One (either pulley or cutting) load acts on 
the overhanging portion of the shaft while 
the other (gear) load acts in between the 
bearings. 

b. Eoth the loads act on the gears and lie in 
between the bearings* 

c* Both the loads act on the gears but one of 
them lies in the over-hanging portion. 

The above three loading conditions have been taken 
into considerations in the Computer programme for bending 
moment and deflection calculations. 

The shafts are designed first for adequate bending 
stress and then checked for deflection. The deflection at 
the gear point is not allowed to exceed a value of m/50, 

B ,2 Bearings: 

The recommended working life of a bearing is about 
eight to ten thousand working hours. There may be several 
types of bearings available for this life and for the condi- 
tions of loading which develop maximum thrust on the bearipgs* 
Therefore, it is recommended here that factors like static 
and dynamic capacity of bearings be tabulated. This helps 
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the designers to select the best type from the available 
stock. 


The relationship between the bearing life and the 
basic dynamic capacity is given as 

L = (C/P) P (b.D 

where 

L - Nominal life in millions of revolutions, 

C - Basic Dynamic Capacity 

P - equivalent radial load in the bearing and 
p - 3 or 10/3 for ball or roller bearings respec- 
tively. 

The relation between life in millions of revolutions L and 
life in working hours is expressed by 

60 * n • Lv, 

II / -n O \ 


where n is the r.p.m. 


For a life of 10,000 working hours 
L — 0 * 6 *n 

and the dynamic capacity for the ball bearing is 

* r 1/p 
C as P . L 

0.333 

= p ( 0.6. n ) 

The equivalent load P is not the maximum load amon- 
gst all the conditions, of loading of the shaft because the 
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gear box is not expected to run always at that speed of the 
spindle. P is function of the individual loads in each type 
of loading, corresponding shaft speeds and the percentage 
time for which they act. In absence of informations about 
the time -sharing of the spindle speeds, average of all the 
dynamic capacities obtained in different loading conditions 
is tabulated here in order to avoid overdesigned bearings. 

B.3 Clutches s 

Three types of clutches are generally used in the 
gear boxes. They are : 

a. Mechanical or Jaw Clutches, 

b. Hydraulic Clutch and 

c. Electromagnetic Clutch. 

The dimensions of the clutches, in majority of designs, are 
not in exact conformity with the requirements of the gear 
box but the clutches are only selected from the available 
stock. Therefore , unless the thermal utilisation factor of 
the selected clutches is more than 10C$, the selection is 
approved. Care should be exercised to select a clutch not 
very much overdesigned, because it costs more and occupies 
more space. Only electromagnetic clutches are discussed in 

the present work. 

Fig. B3 shows an electromagnetic clutch. Either 
the magnetic body or the armature is coupled to the gear of 
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;he drive and rotates on the shaft freely. The remaining 
jart is keyed to the shaft. The preliminary selection of 
;he clutches are done on the hasis of maximum torque to he 
;ransmltted from a gear. 

3.3.1 Thermal Utilisation of the Clutch : 


This is the ratio of the maximum energy absorbed 
oy the clutch to the equivalent capacity of the clutch or 
Th.U . = E/A 0 

who re 

E * maximum energy absorbed and 

A = equivalent capacity. 



where 



equivalent rotational ineitia of the 
system about any shaft , 


y s final angular velocity and 

a 


Q X 


Initial angular velocity. 


To calculate maximum energy 
ray diagram is considered to 


of the clutch, that path in the 
r which different intermediate 


shafts and spindles accelerate to maximum speed from rest 

hen clutch is energised. One such path of the ray diagram 

for the first clutch in Figl. 3 is shown in Fig. 1.4 by 
thick line. 

B.3,2 Angular Velocity of Different Components: 


All the components on a shaft do not rotate at the 
same angular velocity. This is illustrated by Fig. B.4. 
k and A represent the magnetic body and armature respectively. 
The magnetic body of the clutches 1 and 2 and the gears fixed 
to them are free to rotate on the shafts. When first clutch 


is engaged, the components S n , m , A , and rotate with 
speed with the shaft. But the idling parts Mg and 
rotate at a speed of S where 


S 



U 

u 


11 

12 



e 


Therefore total 
from rest is 


energy required to accelerate the 1st shaft 

* < > 

£ (I < + I g e 2 . U^) (B »4) 

( h + % • * 2 } 

2 

^ ^ 11 ^ 1 ^ e qui valient 
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2 o 

e is the factor by which the (GD ) values of idliy rota- 
ting part should be multiplied to give same speed of all 
the components on that shaft. 

Such considerations have been taken into account 
in programming. 

B.4 Design of Keys ; 

Keys are used to fix different components to the 
shaft. In large number of cases the predominant cause of 
failing of keys is crushing. 

The crushing area Ab is given by 
F r p d 

Ab = (B.5) 

d . c 

where 

F = Maximum gear load 

r 

p d = pitch circle dia. of the gear 

d = dia. of the shaft. 

= allowable stress in crushing, 
c 

Only A b is tabulated here and helps the designer to choose 
length and width of the key depending on the suitable type 
of the key. 

B.5 Correction of Gear Tooth : 

The n umb er of gear teeth which have been obtained 
from the optimisation programme shall have non integer values 
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in general* They must be rounded off to the next integer 
with a view that the final speeds at spindle do not change 
by more than 2 From strength considerations, if the pin- 
ions are given positive corrections, their strength increases 
at the cost of already overdesigned gears «- Total correction 
Xrj, required to accomodate the modified number of teeth is 
given as 


t 



Inv (c< ) 


Inv (<^ 0 ) . ( Z gl , + \ ) 

2 tan (ot' 0 ) 

0 (B .6) 


where 

i 

Z = modified number of teeth in pinion. 

cL 

f 

Z^ = modified number of teeth in gear. 


and 


where 


ok 0 = Initial pressure angle 


C< = modified pressure angle 
* tan (_l-n_£L ) 


Q 


* Zb' > 
Z + Zv 


. tan 0( 


a 


unmodified number of teeth in pinion 
unmodified number of teeth in gear. 


B .6 Inspection Data : 

Depending on number of teeth, pitch diameter and 
quality of the gear the inspection data are generated 

1 ationships. 
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APPENDIX 0 


COMPUTER PROGRAMME 

This appendix gives the descriptions of fourteen 
subroutines which add up to make the optimised design of 
gear box complete. A listing of each subroutine and a 
general flow diagram is also included. The input and output 
formats are described and a part of output listing is printed. 

C.l The Programme: 

The programming has been split in two phases. 
Optimisation and Design. The first part is carried out in 
main routine which requires another seven subroutines to reach 
its end. Second phase is carried out in the subroutine "Design 11 
which needs the optimum design vector from Main as input in- 
formation and requires yet another five subroutines for 
component design. 

The purpose of different subroutines together with 
the necessary input for them follows in the subsequent para- 

graphs . 

C.1.1 Main *• 

Purposes To read the input data as given in Art. A.4.2 
To extrapolate the Design vector for a new 
starting point (to be used in subroutine 
"OTCOHS") and to stop iterations when conver- 
gence criteria are satisfied. 



C.1.2 Name of the Subroutine RFCG 


10T 

Purpose: To recognise the groups which can have inde- 
pendent number of teeth and module value for 
one of the driver gears and to calculate the 
factors which are needed in the calculation 
of transmission ratios in case there happens 
to be two composite gear in the first group. 

Input : The configuration of the gear box read in the 
main viz. groups, speeds, composite gears, 
back gears are transferred to this subroutine 
by main through co-mnon statements. 

C.1.3 Name of the subroutine - iiDJST (X) 

Purpose : To correspond different design variables Xj_ s 
in the relevant elements of the matrices 
representing transmission ratios, number of 
teeth and module values and to remodify diffe 
rent elements of the matrices from the compa- 
tibility considerations. 

This subroutine is called every time when 
the change in design vector takes place. 

Input : First input is present design vector. The 
configuration of the gear box transferred 
by the common statement and output from rou- 
tine "KFCG" are also the input informations. 
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C.1.4 Name of the subroutine FSNAL (FO, FT, FP, X,N, R, CSC) 

Purpose : To calculate the value of the objective func- 
tion and to evaluate the normalised constra- 
ints as specified in the article (C.2), To 

r> /*\ 

set the value of penalty function as ICT 
"whenever a constraint is violated (this iden- 
tifies a acceptable solution from unaccept- 
able one). 

Input ; The present design vector, number of variables 
and the value of r in the unconstrained mini- 
misation are the inputs to this subroutine. 
For every evaluation of constraint the sub- 
routine first calls the subroutine '‘ADJST'* 
to modify different matrices, used in the 
constraint evaluatidn. 

C.1.5 Name of subroutine UNCONS (X, N, R) 

Purpose : To carryout unconstrained minimisation for 
any r value in the sequence of r minimisa- 
tion, In process, to identify a situation 
as whether to modify H-matrix or not and to 
stop iteration when the convergence criteria 
are satisfied. 



Input : This subroutine is called from Main every- 
time when r is reduced. Starting design 
vector, number of design variables and the 
value of r are the input informations. 

This subroutine uses three other subrout- 
ines, namely PENAL, SEARCH, HaftlBI at 
different stages. Convergence limits are 
transfered from main through common 
statements. 

,1.6 Name of the subroutine, SEARCH (SIG, X, R, N) 

Purpose : To carry out one dimensional minimisation 
of P (<T ) and evaluate a step size which 
makes the penalty function minimum in the 
present direction of move , using third order 
interpolation. This subroutine is called 
from Main each time a new direction is found 

out. 

Input : The new direction and the present value of 
design vector are input here. Convergence 
criterion has been permanently defined in 
the subroutine itself. 

C ,1.7 Name of the Subroutine GRAD (G, X, N, R) 

purpose : To calculate the gradient at any point m 
the design space by forward interpolation. 



This is called from UNCONS every time a 
step in the old direction is taken and a 
new Design vector is reached. , 

Input : Inputs are the present value of X and number 
of design variables. The step length used 
in the forward difference scheme is perma- 
nently defined in the subroutine ■> 

C ,'ioS Name of the subroutine HmllRI (H, SIG, N) 

Purpose : To modify the Hessian matrix used to calcu- 
late the direction of move at every point. 
Every time step size in the previous direc- 
tion meets the convergence requirements, 
this subroutine is called from UNCONS. 

Input 5 Inputs for the subroutine are the present 
Hessian matrix, the value of step in the 
previous move, the difference of the present 
and previous gradient vector and the number 
of variables. 

0,1.9 Name of the Subroutine : DESIGN 

Purpose : To round-off the modules of gears in diffe- 
rent groups to the next standard module 
available. Round off the number of teeth 
of different gears to make it an integer. 
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Distribute corrections on the gears and 
generate inspection data . To calculate 
percentage deviation of each output speed 
from the desired speed and generate design 
data for clutches, shafts, bearings, keys. 

Input i The final matrices of number of teeth, mo- 
dules and transmission ratios obtained at 
the optimum design variable is input in 
the programme. 

C.1,10 Name of subroutine : CORKBC (ZP,ZB , TRANS ,XT ,XP,XG) 


Purpose : To distribute the total correction on the 
penion and gear for achieving the compati- 
bility. The informations about the graphs 
in I.S. (3756 ) used for this purpose is 
permanently defined in the subroutine in 
Data statement. This subroutine is called 
for correction of each pair of gears* 


Input 


: The total correction value obtained from 
Design, rounded number of teeth in gears 
and pinions and the transmission ratio 
between them. 


C.1,11 Name of the subroutine - INSFSC (Z,DP,X,V,Am,IQLTY) 


Purpose : To calculate, the number of teeth over which 



the measurement is to "be taken and to print, 
upper, lower and nominal value of tolerences 
for each gear. 


Input 

C.1.12 Name 
Purpose 


Input 


: The input informations are the number of 
gear, its diameter. Corrections on it and 
the quality of the gear. 

of subroutine : SHBG 

: To obtain shaft diameters and basic dynamic 
capacities of all the bearings used in the 
gear box considering all possible loading 
arrangements of the shaft. To calculate 
max imum bending moment and deflection under 
the worst condition of loading for each 
shaft . 

: The input informations are the approximate 
distance of gears from left hand bearing, 
the number of driver gears in each group, 
belt and cutting forces, Young 
elasticity, allowable value of stresses 
(in bending and torsion), efficiency of each 
drive, approximate lengths of the shafts 
and the life of bearings desired. 



C.1.13 Name of subroutine : SINGLE (Fl, F2, XL, LISl, DIS2, 

AN, TORQ, ALH) 

Purpose : To calculate the equivalent bending moment, 
deflection and the basic dynamic capacities 
of the bearing in a case when shaft is loa- 
ded by an overhung force. 

Input : Magnitudes and points of application of the 
driver and the driven forces, angle between 
the forces, torque required to be transmitted 
and the working life of bearing. 

C.1.14 Name of subroutine : CLUTCH 

Purpose : To select a suitable clutch at the required 
places and to calculate the thermal utilisa- 
tion of the clutch selected. 

Input i Input for the subroutine is the number of 
teeth, module values and face width of the 
gears, length of shafts, specific weight of 
the material, speed diagram, indication of 
the portion (magnetic body or armature) 
attached to a gear and manufacturer's data 
for the clutches. These are permanently 
read in the clutch programme in data state- 


ments 



C .2 Input to the main programme : 


The required input data is described below in the 
order in which it is read. The Fortran symbolic representa- 
tion is given on the left and an explanation on the right. 


Identifier 

NSPD 

NSH 

IBGSH 

ILAP 


NGMX 


IQLTI 

HP 

ETA 


Explanation 

_ Total number of spindle speeds. 

- Number of shafts in the gear box. 

- Name of back gear shaft. If there 
is no back gear shaft then IBGSH=0. 

_ Number of steps by which the lowest 
speed of a shaft preceding back 
gear shaft is lagging to the speed 
transfered to the same shaft by 
clutch coupling. (This is 4 in 
Fig. 1.4) . Put zero if no back 
gear shaft is present. 

- Maximum number of gears that is 
possible on any shaft. Count 2 in 
case there is a composite gear. 

- Quality of the gears desired to be 
used in the main drive. 

Hated horse power capacity. 

- Efficiency of the drive at each 


stage. 
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Identifier 

Explanation 

RPMI 

- Input Motor speed. 

REML 

- Lowest desired spindle speed. 

RPMH 

- Highest desired spindle speed. 

TMAX 

- Maximum allowable torque. 

IP ( ) 

- Number of pairs of gears in each 

group . 

ICAR ( ) 

- Characteristic of each group. 

ICOMP ( , ) 

- Composite gears 

= 0 if the driver er the driven 

gear is not a composite gear. 

= I where I is the order number 

of the transmission it produces. 

.NOCL . 

_ - s .TRUE, if no clutch in the drive 

= .FALSE, if there is clutch. 

ICL ( , ) 

- Read only when NOCL = .FALSE. 

= Q when gear is not attached to 

the clutch. 

s i if armature is keyed to the 

gear. 

s 2 if magnetic body is keyed to 

this gear. 

ANG ( ) 

^ Included angle between two adjacent 

shafts . 

ALSH ( ) 

- Approximate length of the shafts. 
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Identifier 


Explanation 

DIS ( , ) 

- 

Approximate distances of the gears 

from left hand bearing. 

GDCH 

- 

p 

GD - value of chuck. 

GDP 

- 

GD 2 - value of Pulley. 

RHO 

- 

Specific weight of carbon steel. 

PLF 

- 

Pulley force. 

CUTF 

- 

Maximum cutting force. 

DSP 

— 

Distance of the Pulley from left 

hand bearing. 

DSJ 

— 

Distance of the job from right hand 

bearing. 

AAMOD ( ) 

— 

A column of standard module-value 

from 1 to 10 a 

SC 

- 

Material factor 

HOUR 


Hours per day each gear is expected 

to remain active. (Approximate 

value ) . 

CF 

*• 

Factor by which the values of r 

should be divided each time. 

SENS 

• 

Sensitivity factor for the objective 

function. 

N 

« 

Number of independent design 


variables. 
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Identifier 


Explanation 


N 

XmlN ( ) 

SPAN ( ) 

XI ( ) 

SIGT 

SIGC 

si'IGB 

SiGtf 

FSC 

FNUT 

FS 

E 

ALH 

EPSl 

EPSS 

EPS3 

EPS4 

ITIR 


- Number of independent design variables - 

- Minimum value of independent design 
variables . 

- Maximum spread of independent des- 
ign variables. 

- Independent design variable (initial 
value ) . 

- Allowable stress in torsion. 

- Allowable stress in crushing. 

- Allowable stress in bending. 

- Allowable working stress. 

- Factor, stress concentration. 

- Nutsc. factor 

- Factor of safety. 

- Young's modulus 

- Desired bearing life in hours. 

- Required accuracy in penalty functio.j. 

- Required accuracy in design 
variables. 

- Required accuracy in one dimen- 
sional minimisation. 

- Required accuracy in the object 
function. 

- Number of iterations in -minimisa- 
tion for a particular r. 
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C.3 Programme Output: 

The input data is written out in the same format 
in which it is read. 

Other outputs are presented with self explaining 

headings. 

C .4 Programme listing and sample output; 

The subsequent pages contain the listing of all 
the subroutines and a portion of output. 

























CXj Cvl 
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i J 0 B 

'FTC MAIN 

LOGICAL NOCL 

C0MM0N/LAB2/NSPD*NSH,IBGSH*ILAP,N6MAX 

C0MM0N/LA63/HP»ETA»RPMI 9 RPML » RPMH , TMAX ,PH I 

COMMON /LAB5/ IP, ZS,SPD 

COMMON/LABA/ICL, ICAR 

COMMON /LAE'?/ AM »TR ,DSH, ALSH 

COMMON/LABIO/ANG.DIS 

COMMON/LAB13/ ICOMPsSPDM 

COMMON/LAB15/ CQRR 

COMMON /L AB32/SP AN, XM IN ,SENS 

C0MM0N/LAB34/NC 

COMMON/LAB35/CZ, SC, RTL 

COMMON /LAB42/G 

C0MM0N/LAB44/CEN ,AREA 

COMMON /LAB46/AAM00 

C0MM0N/LAB47/ EPPS 1 ,EPS2 ,EPS3 , EPS 4 , I T I R 
COMMON /LAB 5 C 7 SIG.T ,S I GB ,S IGC 

DIMENS ION AM ( 6 , 6 ) » T*? { 6 ,6 ) ,CEN ( 6 ) ,AAMOD I 24 ) , SPD ( 6 , 6 ) ,S^DM ( 6 , 6 ), 

1 IP (6) , DSH ( 6 ) , I COMP (6,6) , ZS ( 6 , 6 ) , AR EA { 6 , 6 ) , ALSH (6) ,AMG(6) »DI S ( 6 ,6 ) , I CL ( 

2 I CL (3, 6) , ICAR ( 6 ) 9 CORRC 656 ) 

DIMENSION X ( 1 1 ) , X 1 ( 1 1 > , XR ( 1 1 ) * XL ( 1 1 ) ,X0(11) , GX(100),G(11) 

DTMFNSION XMIN (11) 9 SPAN (11 ) ,XM( 11) 

6 format (/, i x,* extrapolated values of variables at next r-value*> 

l/, 2 {i X,5F17 o6/) ) 

17 FORMAT ( IX , 616 ) 

3 FORMAT (///,1 X,* NORMALISED OPTIMUM DESIGN VARIABLES*,/ 

12(1 X 9 5 F 1 C 6 / ) ) 

5 FORMAT ( / * 1 X,* ACTUAL DESIGN VARI ABLES* 9 / , 2 ( 1 X 95 F 1 O 04 /)) 

27 F0RMAT(1X96F10o2) 

37 F0RMAT(1X»3F12.6) 

•4 FORMAT (* INITIAL VALUE OF R=*,F10*2) 

77 FORMAT (1X.6F11.3) 

87 F0RMAT(1X94F9 o 2) 

97 FORMAT (IX 96 F 80 I) 

117 F0RMAT(lX,10F6ol) 

L 27 FORMAT (12F6o3) 

137 FORMAT ( IX 96 I 6 ) & 

147 FORMAT ( IX » 15 14) 

157 FORMAT (1X93012) 

167 FORMAT (IX, 13) 

177 FORMAT (lX,3F5o2) 

187 FORMAT (1X,9F6.2) 

197 FORMAT (IX, LI) 

07 FORMAT (1X,2F8o3) 

17 FORMAT (1X,4F7 0 3,I4) 

227 FORMAT (1X,2F5*1) 

500 FORMAT (* OBJECT FUNCTION *9F10„2,/ 

1* PENALTY * , F 1 0 e 2 , / , 

2* TOTAL FUNCTION*, F10o2) 

501 FORMAT (1 X,* DESIGN VAR I ABLES* , / , 2 ( l X»5FlOo6/)> 

502 FORMAT (1 X,* CONSTRAI NTS* , / , 20 ( 1 X»5FlOo6/)) 

504 FORMAT f IX *80 ( 1H— > ) 

503 FORMAT (1 X,* GRADIENT VECTOR* 9 / , 2 ( 1 X»5F10«6/)) 

506 FORMAT (1X»*NUMBER OF TEETH IN GEARS* » /IX, ( 1X,6F10*2/ > > 
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507 FORMAT (1X»*L0EST TRANSMISSIONS FROM THE GROUPS* * / 1X»6F10»4) 

508 FORMAT (IX,* CENTRE DISTANCE**, 6F8 <>2 ) 

READ 17,NSPD*NSH» I BGSH , I LAP , NGMAX , I QLTY 
PRINT17»NSPD»NSH» IBGSH, ILAP*NGMAX,IQLTY 
READ 27 ,HP ,ETA,RPMI , RPML , R PMH » TMAX 

PPI NT27.HP, ETA, RPMI , RPML ,RPMH » TMAX 
READ 137, ( T P ( I ) ,1 = 1, NSH) 

PR INTI 37 , ( IP ( I) , 1=1, NSH) 

READ 137 , (ICAP( I ) ,1 = 1, NSH) 

PRINT 137. ( ICAR ( I ) ,1=1, NSH) 

READ 157, U I COMP ( I , J) , J = 1 , NGMAX ) » 1=1 , NSH ) 

PR I NT 157 , ( ( I COMP ( I , J) » J = 1 , NGMAX) , 1 = 1 , NSH) 

READ 1 97 , NOCL 
PR I NT 1 97 , NOCL 
IF (NOCL) GO TO 237 

READ 147, U ICL( I *J> ,J = 1, NGMAX) ,1=1,3) 

PR I NT 147 , e ( I CL ( I > J ) »J = 1> NGMAX ) ,1 = 1,3) 

237 CONTINUE 

READ 97, (ANG( I ) , 1=1 ,NSH) 

PRINT97, (ANG( I ) , 1=1, NSH) 

READ 97, (ALSHt I ) ,1=1, NSH) 

PRINT97, <A(_SH( I ) ,1 = 1, NSH) 

READ 117, ( (DIS( I , J ) , J = 1 , NGMAX ) , 1 = 1 ,NSH ) 

PRINT117»( ( D I S ( I , J ) ,J=1 , NGMAX ), 1=1, NSH) 

RE AD 37 » GDCH , GDP , RHO 
PR I NT37 , GDCH ,GDP , RHO 
PEAD 87,PLF,CUTF,DSP,DSJ 
PPINTS7,?LF ,CUTF,DSP,DSJ 
READ 127, (AAMOD( I ) ,1 = 1,24) 

PRINT 127, ( AAMODt I ) ,1=1,24) 

READ 207, SC, HOUR 
PRINT207 , SC, HOUR 
PEAD227 , CF .SENS 
PRINT 227, CF .SENS 
READ 167, N 
PRINT 167, N 

READ 187, (XMIN( I ) ,I=1,N) , ( SPAN ( I) ,1 = 1 ,N) 

PR INTI 87, (XMINC I ) ,I = 1,N) , (SPAN(I) , I = 1,N) 

READ 1 87 , ( XI ( I ) ,1*1, N) * 

PRINT187»(X] (I) » I = 1 »N ) 

READ 77,SIGW,FSC,FNUT,FS,E,ALH 

PRINT77»SIGW,FSC,FNUT,FS,E,ALH 

READ 1 77,SIGT,SIGC,SIGB 

PRINT1 77,SIGT,SIGC,SIGB 

READ 217,EPS1,EPS2,EPS3,EPSA, ITIR 

PRINT217,E d S1,EPS2,EPS3,EPS4,ITIR 

RTL = RPML/RPMI 

MGRP = NSH-1 

AMSPD = NSPD 

ANSPD = AN$PD»1 « 0 

T = 1 o 0/ ANSPD 

PHI = (RPMH/RPMLT**T 

c 2 = 0 f> 3442*( H 0ui®b o 821)**0 o 218-l o 385 
CALL RECG 
DO 505 I = 1 ,N 

505 X l ( I ) = ( XI ( I )-XMIN ( I) ) /SPAN ( I ) 

CALL PENAL ( FO , FP , FT , Xl »N , 1 « 0 ,GX ) 
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IF (FT oGT ol0o0E+20) GO TO 448 
DO 43 T = l , NC 

43 IF(GX(I).LT.0o 2) FP = FP+5 a 0—1 o 0/GX ( I ) 

PRINT 500 , FO ? FP ? FT 
PRINT 502, (GX( I ) , I=1,NC) 

R = FO/FP 
PRINT 44 , P 
L3 = 0 

I CALL UNCON. 5 ( XI »N » R ) 

DO 5 I = 1 ,M 

5 XN ( I ) = XI ( I ) 

CALL PENAL ( FO , FP , FT ,XM *N , R ,6X ) * 

PRINT 500 , FO , FP , FT 
PRINT 501 , ( XM( I ) » 1=1, N) 

PRINT 502, ( G X ( I ) , I = 1,NC) 

PRINT 503 , (G( I ) , 1=1 »N ) 

PRINT *506 »( ( Z 5 ( I » J ) , J = 1 , NGMAX ) ,I*1»NSH) 
PRINT 507 , ( TR ( I , 1 ) , 1 = 1 »MGRP ) 

PRINT 50 8 , (GEN ( I ) 9 1 = 1 , N6RP ) 

PPIMT 504 
L 3 = L3 + 1 

IF (L3.EQ.1) GO TO 6 

I F ( ABSt (FC-FLO) /FO) oLT„EPS4) GO TO 50 
O SS = ABS( (FT-FL) /FL) 

I F ( SSoGT o E D S1 ) GO TO 3 
DO 10 I = 1 , N 

A 1 * ABSt (XMm-XL(I) )/XL(I) ) 

IF( Al.GToEPS2» GO TO 3 
10 CONTINUE 

DO 11 I = 1 , N 

II X0(I) = (CF**0o5*XM(I)-XL(I))/(CF**0. 5-1.0) 
CALL P ENAL ( FO , FP , FT 1 , XO ,N , R , 6X ) 

IF (FTloGT .10o0E+20) GO TO 3 
IF (FTloGT oFT) GO TO 3 
SS = ABS ( (FTl-FT) /FT) 

IF (SSoGToOoOOl) GO TO 3 
DO 12 I = 1 , N 

A 1 = ABS ( ( XO ( I ) — XM ( I ) ) /XM ( I ) ) 

IF (A1.GT.0.C01) GO TO 3 
12 CONTINUE 

GO TO 26 

6 RL = R 

DO 7 I = 1 , N 

7 XL ( I ) * XN( I > 

FLO * FO 

FL = FT 
R = R/CF 
DO 8 I = 1 , N 

8 XI ( I ) = XN( I > 

GO TO 1 

3 RL = R 

FL = FT 
FLO = FO 
R = R/CF 
DO 15 I = 


1 , N 
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XP(I) = XM( I )+l„ 0 /CF** 0 . 5 *(XM{ I )— XL( I ) ) 
PRT^T 16 , (XR( I ) , 1 = 1 , N) 

C;Li. PFNAL (F 0 ,FP,FT 2 ,XR,N,R,GX) 

It- ( FT 2 .GT o lCoOE +20 ) GO TO 18 

no ] o j = i , is; 

XI ( I ) = xp ( I ) 

GO TO 1 
DO 21 I = 1 ,N 
XI ( I) = XM ( T ) 

GO TO 1 
CONTINUE 
DO 51 I = 1 * N 
XO ( I ) = XI ( I ) 

PRINT 23 »(XO(I)»I=l»N) 

DO 24 I = l,N 

XO(I) = XCtI>*SPAN(I}+XMlN(I) 

PRINT 25 , (XO<I ) , 1 = 1 ,N) 

CALL DESIGN 

STOP 

END 


FTC ADjST 

SUBPOUTINE AD JS T I X ) 

DIMENSION ID1 (2 ) , ID2(2 ) ,IK(6) ,X( 11) 

DIMENSION AM ( 6 , 6 ) ,TR(6,6) ,CEN(6) , AAMOD ( 2.4 ) , SPD ( 6 ,6 ) , SPDM ( 6 ,6 ) , 

1IP( 6) , DSH ( 6 ) , I COMP (6,6 ) , ZS C 6 , 6 ) , AREA ( 6 ,6 ) »ALSH(6> , ANG f 6 ) ,DIS(6»6),ICL(0 ? 
2 I CL ( 3 , 6 ) , ICAR (6) ,C0RR(6,6) 

COMMON /LAB2 / MSPD » NSH , I BGSH , ILAP ,NGMAX 

COMMON/LAB VHP »ETA,RPM I > PPML , RPMH , TMAX , PH I 

COMMON /LAE 6/ I P,ZS, SPD 

C0MM0N/LAB6/ICL, ICAR 

COMMON/LAB9/AM,TR,DSH,ALSH 

COMMON/LABIO/ANGjDIS 

COMMON/LAB13/ I COMP, SPDM 

COMMON /LABI 5 / CORR 

COMMON /LAB 17 /MX 1 »MX2»FS1 > I CTWO , I TV , I K , I Dl , ID2 
COMMON/LAB34/NC 
COMMON /LAB3 5 /CZ ,SC,RTl 
COMMON /LAB33/FO 
COMMON /LAB44/CEN, AREA 
C0MM0N/LAB46/ AAMOD 
NGRP = NSH-1 
DO 46 I = 1 , NSH 
DO 46 J = 1 » NGMAX 
ZS(I,J) ' Ool 
T P ( I , J ) = 0.1 
AM ( I * J ) = 0.1 
SPD ( I » J ) =0.1 
S P D M ( I » J ) * 0.1 
AREA( I , J ) * 0.1 
46 CONTINUE 

DO 100 'I = 1 , 1 TV 
K = I K < I ) 

ZS ( K ♦ 1 ) = X 1 2*1-1 ) 


t AM ( K » 1 ) = X ( 2*1 ) 

NDV = ITV*2+1 
I F ( ICTW0.EQo0)G0 TO 110 
SI = X ( NDV ) *FS1 

TR ( 1 » 1 ) * ( S 1*PH I **MX2* ( PH I**MX1-1 0 0 ) — ( PH I**MX2— 1 ,0) ) / (Phi**MX 2-PH!**MX' 
II **MX1 ) /PHT**(ID1 ( 1 )-l ) 

TR ( 2 , 1 ) = SI /TR ( 1 » 1 ) 

GO TO 120 

0 TR ( 1 » 1 ) = X(NDV) 

NDV = NDV+1 
TP ( 2 » 1 ) = X ( NDV ) 

0 DO 130 i = 3 , NGRP 

' IF(I.FOoNGRP) GO TO 140 
IF { I c, EQcTBGSH) GO TO 15 0 
. ' NDV = NDV+1 

TP ( 1,1) = X ( NDV ) 

GO TO 130 

T R ( I » 1 ) ~ 1 » 0 /TR ( I — 1 , 1 ) /PHI ** I LAP 
O0 TO 130 
m 3 SUV = loO 

NGRP1 = NGRP-1 
DO 160 J = 1 » MGR PI 
>0 SUM = SUM*TR(J.l) 

TR ( 1,1) = RTL/SUM 
30 CONTINUE 
BB = RPM I 
AA = RPM I 

DSH(l) = l„25*(716o0*HP*1000 o 0/AA/2o >**0.3333 

DO 169 I = 1 , MGRP 
K = IP ( I ) 

DO 17C J = 1 , K 
K 1 = IP(I+1)+J 

TP(I,J) = TR £ I ,1 )*PHI**( ICAR( I )*(J-1) ) 

* Tp ( 1 + 1 ,K1 ) = loO/TR ( I » J ) 

SPD ( I , J) = AA 

SPDM( I , J > = BB 

SPD ( I + l.Kl) = A A*TR ( I ,J> 

L 7 0 SPDM(I+1,K1) = BB*TR(I,J) 

AA = AA*TR(I,1) 

TORO = 716,0*HP*1000»0/AA*FTA**I 
IF { TORQoGT oTMAX ) TORO * TMAX 
DSH ( I + 1 > = 1.25*(TORQ/2. )**0o3333 

169 BB = BB*TP ( I , K ) 

DO 180 I * 1 , NGRP 
K = I P t I ) 

IF (IoNEoIBGSH) GO TO 220 
CEN( I ) = CEN( I — 1 > 

AM ( I , 1 ) = AM ( 1-1,1) 

DO 230 J = 1»K 
J1 = IPU+U+J 

ZS(I»J> = 2oO*CEN(I)*TR(I,J)/AM(I,l)/(TP(I,J)+IoO> 

2S ( 1+1 » J1 ) = ZS( I ,J)/TR( I,J) 

AM ( I , J) = AM ( I ,1 ) 

AM ( I+1,J1 ) * AM ( I , 1 ) 

230 CONTINUE 
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GO TO 180 

11 * 0 

12 = 0 

DO 190 J * 1 »K 
K 1 = IP{I+1)+J 

IF ( I COMP < I >j) 0 NE 0 O) II = I COMP ( I , J ) 

I F ( I COMP ( I +1 »K1 ) o EQ sO ) GO TO 190 
JP = K 1 

12 * I COMP ( I + l.Kl) 

9 CONTINUE 

I F ( TloNEoP) GO TO 180 
21 = ZS(I,1) 

<2 = IP ( 1+11 + 1 

ZS ( 1 + 1 ,K2 ) = ZS(I,l)/TR(I.l) 

Z 2 = ZS ( I +1 » K2 ) 

CEN { I ) = (Z1+Z2 ) * AM ( I , 1 ) /2o0 
DO 200 J = 1,K 
k i = ipn+u+j 

ZS C I » J J = <Zl+22)*Tp(I,J)/(loO+TR(I,J) ) 

ZS { 1+1 »K1 ) = ZS(I,J)/TR(I,J) 

AM ( I s J ) = AM ( 1,1) 

AM ( 1 + 1 ,K1 ) = AMI 1,1) 

10 CONTINUE 

I F t I 2 o EO o 0 ) GO TO 180 
ZS( 1+1,12) = ZS< I+1,JF) 

13 = I P ( 1+2 ) +1 2 

ZS ( 1+2,13) = ZS( I+1,I2)/TR( 1+1,12) 

21 = ZS( 1+1, 12) 

Z 2 = ZS( 1 + 2, 13) 

C E N ( I + 1 ) = (Zl+Z2)*AM(I*l)/2.0 
!< = IP ( 1 + 1 ) 

DO 210 j = 1 , K 
J1 = I P( 1+2 ) +J 

ZS(I + 1,J> = (Z1+Z2 )*TR( I + 1,J) / (1 oO+TR I 1 + 1 , J ) ) 

ZS ( I +2 » J1 ) * ZS( 1 + 1 ,J)/TR{ 1 + 1 ,J) 

AM ( 1+1 , J ) = AM (1,1) 

AM ( 1 + 2, J1 ) = AM ( I ,1) 

210 CONTINUE 
180 CONTINUE 
FQ = 0 o 0 

DO 45 I = 1 , NGRP 
45 FQ - FQ+CENC I ) 

RETURN 

END 

EIBFTC RECG . 

. • SUBROUTINE. RECG 

£ ' 01 PENSION I K ( 6 ) » I D 1 ( 2 ) »• I •" 2 ( i ) 

DIMENSION X*l6>6) , T P ( 6 * 6 ) >CEN<6> »AA"'0D<24) >.Spp(6,s) >Sdqm( 6 , 6) , 

U P « 6 ) ' DSh ( 6 ) »I C0M p ( 6*6 ) »ZS ( 6 , 6 ) » AREA. ( 6 , 6 ) , ALSh f 6 ) ♦ A.2G?6 ) » D I S ( $ , 6 ) » I CL ( 3 , 
2 I CL ( 3 i- 6 ) ? I CAR ( 6 ) > CORR ( 6 > 6 ) 

COMMON /LAE2/NSPD»NSH, I BGSH , I LAP , NGMAX 
C0MM0N/LAE3/HP»ETA»RPMI ,rpml»rpmh,tmax,phi 
COMMON /LAE 5/ IP»ZS,SPD 
C0MM0N/LAB6/ ICL , I CAR 
COMMON /LAB 9 /AM , TR »DSH , ALSH 



COMMON/LABl-'VANG>DIS 
C0MM0N/LAB13/ ICOMP,SPDM 
C0MM0N/LAB15/ CORR 

COMMON/LABl7/MXl»MX2.FSlt ICTWO»ITV» IIC*ID1 # ID2 

C0MM0N/LAB44/CEN»AREA 

COMMON /L A B46/AA MOD 

MGRP = NSH-1 

I TV = 0 
ICTWO = 0 
NDV = 0 

30 100 I = 1 9 NGRP 

II - I P ( I ) 

I F ( I o EQ o I BGSH } GO TO 100 
I Nl n G = 0 
DO 110 j = 1,11 

110 IF( ICOMPC I .NEoO) INDG = INDG+1 
IF (INDGoNE.O) GO TO 100 
NDV = NDV +1 

I TV = ITV +1 
IK(NDV) = I 

100 continue 

IF< IP(1) o LT o 2o0R o IP(2) »LT „ 2 ) RETURN 

II a 0 

12 = IP(1) 

DO 120 J = 1,12 
J1 = IP(2)+J 

IF (IC0MP(2»J1).EQ.0) GO TO 120 
II = 1 1+1 
IDl(Il) = J 
I D2 ( 1 1 ) = I COMP ( 2 » J1 ) 

120 CONTINUE 

I c ( 1 1 o NEc 2 ) RETURN 
ICTWO = 1 

MX1 = ICAR(1)*IABS<ID1(1)-ID1(2) ) 

MX2 = ICAR (2)*I ABS( ID2 ( 1 )-ID2<2 > ) 

RETURN 

END 

SIBFTC SUB3 

SUBROUTINE HMATRI (H»SIG»N) 


SIBFTC SEARCH 

SUBROUTINE SEARCHCSIG,X,N,R) 

C0MM0N/LAB34/NC 

COMMON/LAB23/D 

DIMENSION X ( 9 ) » X 1 ( 9 ) »D(9) »GX 1(100) » GX2 ( 100 ) 

COMMON /LAB 24 /GX 1 >GX2 

JJ = 0 

II = 0 

ST = OoO 

DO 1 I * 1»N 

XI ( I ) = X ( I ) 

IF (ABStDC I ) ) oGToST) ST = ABSlDd)) 

1 CONTINUE 

ST = 0.001/ST 

CALL PENAL (F0»FP>FT1»X»N»R»GX1) 
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DO 2 I = 1 » N 
X { I ) = XI { I )+D( I )*ST 

2 CONTINUE 

CALL PENAL ( F01 » FP » FT >X?N>RsGX2) 

PTC = 0 o 0 2 * F 0 
GA = FOl/ST-FO/ST 
A = OoO 
DO 3 I = 1 » NC 

3 GA = GA-< (GX2( I )-GXl ( I ) )/ST) /GX1 ( I)**2 
D2 = -PTC*2.0/GA 

D2 = AMINl(leO»D2) 

5 LI = 0 

DO 20 I = 1 » N 
X ( I > = XI ( I )+D2*D( I ) 

IF (X( I)oLToOoO.OR.X(I) .GTol.O) L1 = L1 + 1 
2 ° CONTINUE 

IF (LloNEoO) GO TO 40 
C - L !.. PENAL (F01»FP,FT2,X,N*R»GX1) 

r = 02 

IF ( ~T2 -LT.lOeOE+20 ) GO TO 18 

; - -■ / 2 o o 

T- ■ Dl.cr^f .000001) GO TO 5 
s I G = 0 -- 0 
PET U P N 

-+ DO ]3 I = 1,M 
X ( I ) = X 1 ( I ) 

13 X ( I ) = XI ( I ) +02#D ( I ) 

CALL PENAL (F01,FP*FT2»X»N*R»GX1) 

IF (FT2 o LT,10,0E+20) GO TO 18 
D2 = P2/ 1 o 5 
IF ( D 2 o GT e C ) GO TO 4 
B = C 
GB= GC 
FT 2 = FT 3 
GO TO 8 
18 CONTINUE 

DO 6 I = 1»N 

6 X ( I ) = X( I 5+D( I )*ST 

CALL PENAL ( F02 , FP » FT , X *N » R » GX2 » 

GB = (F02-F01)/ST 

B = ( X ( 1 )-Xl ( 1) ) /DC 1) 

DO 7I--1 ,NC 

7 GB = G8-((GX2(I)~GX1(I))/ST>/GX1(I>**2 
JJ = JJ+1 

IF (JJoNEcl) GO TO 21 
G C = G B 
FT3 = FT2 
C = B 

21 IF ( GBoGToOoOoORoGB.GToGA ) GO TO 8 
GA = GB 
A = B 
FT1 = FT 2 
D2 = 2oO*D2 
GO TO 4 

8 Z = 3.0*(FT1-FT2)/IB~A)+GA+GB 
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'■/ = (Z**2~GA#GB )**0 o 5 

SIG = _ B~< GR+W-Z) / (GB-GA+2oO*W) *(B-A) 
IF { SIGoGTo Ao AND 0 SIG 0 LT 0 B) GO TO 50 
* SIG = A 
GO TO 51 
- CONTINUE 
II = I 1-4-1 
DO 10 I = 1 » N 
X ( I ) = XI ( I )+SIG*D< I ) 

CALL PENAL { F04,FP ,FT4 ,X ,N *R »GX1 > 

DO 11 I = 1 ♦ N 

X ( I ) = X ( I ) +D ( I )*ST 

CALL PENAL ( F05 , FP , FT , X » N , R , GX2 ) 

GC = (F05-F041/ST 
DO 14 I =1 »NC 

GC = GC-( t (1X2 C I )-GXl( I ) }/ST)/GXl(I )**2 

IF (GCoGToOoOoORoGCoGToGA) GO TO 15 

GA = GC 

A = SIG 

FT 1 = FT 4 

C,r> TO 16 

C-P = GC 

8 = SIG 

FT2 = FT 4 

IF ( I I o L T o 2 ) GO TO 8 
CONTINUE 
DO 97 I = 1 » N 
X ( I ) = XI ( I ) 

RETURN 

END 


3FTC PENAL 

SUBROUTINE PENAL ( FO ,FP , FT ,X ,N ,R ,GX ) 

DIMENSION X(ll) » GX ( 100 ) *SPAN(11) »XMI N( 11 ) »X1 ( 11 ) 

DIMENSION AM(6?6)>TR(6»6) >CEN(6) »AAM0D{24) » SPDI 6 *6 ) ?SpDM { 6 *6 > » 

]IP( 6) s DSH ( 6 ) » I COMP (6 , 6) »ZS ( 6 *6 ) » AREA< 6 *6 > * AL-SH ( 6 ) »ANG(6) *DIS ( 6 »6 ) > I CL ( 
2 I CL (3, 6) » ICAR (6) 9CORR(6>6) 

COMMON /LAE 2 /NSPD»NSH, IBGSH 9 ILAP 9 NGMAX 

CO w MON/LAB3/HP 9 ETA 9 RPMI 9 RPML , RPMH, TMAX , PH I 

COMMON /L AB5 / I P 9 ZS 9 SPD 

COMMON/LAB 6 / ICL 9 ICAR 

COMMON/L AB Q / AM 9 TR 9 DSH 9 ALSH 

COMMOI'/LABl O/ANG 9 DIS 

COMMQN/LAB13/ ICOMP 9 SPDM 

COMMON /LABI 5 / CORR 

C 0 V M 0 N / L A B 3 2 / S P A N 9 X M I N 9 S E N S 

COMMON /LAB 34 /NC 

COMMON/LAB35/CZ 9 SC 9 RTL 

COMMON /LAB33/FQ 

C0MM0N/LAB44/CEN 9 AREA 

COMMON /LAB46/AAM0D 

NgRP = NSH-1 

F 1 ( X 9 C ) = X** (-0ol85)*EXP (— C ) 

F2(X) = 0.0294*X-1.3+0o 000169*X**2 
F3 ( X ) = 0 » 222-0 o 00204*X+0 0 0000115*X*#2 
p4(X) = — 0 o 0046 1 +0o000C653*X+0«0000003*X**2 
NC = 0 
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DO 1 I = 1 , N 

MC = NC+1 
GX(NC) = X ( I ) 

XI ( I ) = X ( I ) 

. X ( I > = X1(I>*S PAN ( I )+XMIN( I ) 

CALL ADJST(X) 

DO 2 I = 1 »NGRP 
II = I P ( I ) 

MC = NC+1 

GX(NC) = 1 oO-Oo 25/TP. (1*1) 

MC = NC+1 

GX(NC) = 1„0-TR( I.ID/2.0 
NC = NC+1 

GX(NC) = 1.0-lo75/AM( 1,1) 

MC = NC+1 

GX(NC) = 1 ,0~AM{ I ,1 ) /10e0 
D^ 3 J = 1,11 
J1 = IP(I+1)+J 
MC = NC+1 

GX(NC) = 1 i 0— 1 4 0 0 / Z S ( I+I 9 JI) 

Hr = mc+1 

GXfNC) = 1.0-Z5(I,J)/120.0 
MC = NC+1 

GX(MC) = 1 oO“ZS ( 1+1 , J1 ) /12n„0 
MC = NC+1 

GX(NC) - 3 o 0“DSH( I ) /<0 o 8*(ZS( I ,J)-2o35 )*AM( I »J) ) 

MC = NC+1 

GX(NC) = 1.0*-DSH(I + l)/(0o8*tZS(I+l,JlW.35)*AM(I+l,Jl)) 

NC = NC+1 

GX(MC) = 1 .0-14 0 O/ZS ( I , J > 

YZ = F2(ZS( I + 1,J1) )+F3(ZS(I + l,JD )*ZS( I , J )+F4( ZS ( 1+1 , Jl ) )*ZS(I,J>**2 
1*2 

STK = FI (SPD( I > J) »CZ)*YZ*SC*9.0*AM< I »J)**2.8*SPD( I *J)*ZS(I»J) 

1/25 o 4**0 t>8/lo0E+8*1774 o 0 

NC = NL+1 

GX(NC) = 1 » P— HP /STK 

STK - FI ( S D D ( 1 + 1 , Jl ) ,CZ) *YZ*SC*QoO*AM( 1 + 1 , Jl ) **2 » 8*SPD ( 1 + 1 , Jl )*ZS (1 + 1 
1 1+1 , Jl ) /25,4**0o8/lo0E+8*1774 e n 
NC = NC+1 

GX(NC) = 1,0-HP/STK 

CONTINUE 

CONTINUE 

3 FO r ’N AT ( IX » 5 E15 » 3 ) 

XX = G X ( 1 ) 

DO 8 I = 1 » N C 
I F ( XX 0 GF 0 C-X ( I ) ) XX = GX( I ) 

CONTINUE 

IF(XXoLToO„0) GO TO 9 
SUM = OoO 
DO 10 I = UNC 
0 SUM = SUM+R / G X ( I ) 

FP = SUM 
FO = FQ 
FO= FO**SENS 


5 U 
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~ T = F 0+ F P 
GO TO 11 
FT = in, of +30 
1 DO 3 I = 1 , m 

X ( I ) = X 1 ( I ) 
RETURN 
END 


BFTC GRAD 

SUBROUTINE GRAD(X,N,G,R) 

DIMENSION X ( 1 1 ) »GX 1 ( 100 ) ,GX2 ( 100 ) ,G ( 1 1 ) 

C0MM0N/LAB34/NC 

COMMON /LAB24/GX 1 »GX2 

CALL PENAL(F0,FP,FT,X,N,R,GX1) 

DO 2 I = 1,N 
A = X ( I ) 

ST = AM INl(0o001*A/l000o0) 

X ( I ) = X ( I ) + ST 

CALL PENAL ( F01 , FP , FT , X »N »R ,GX2 ) 

G(I) = ( FC1-F0 ) /ST 
DO 3 J = l, NC 

G ( I ) = G< I )-( (GX2 ( JJ-GX1 ( J) )/ST> /GXK J)**2 
3 CONTINUE 

X ( I ) = X ( I )— ST 
2 CONTINUE 

RETURN 
FND 


EIBFTC UNCONS 

SUBROUTINE UNCONS ( X ,N , R ) 

DIMENSION Y(11)*D{11),GX(100)»X(11)*G(11)»G1<11)*H(11,11» 

COMMON /LA.B3 8 /Y 
COMMON/LAB23/D 
COMMON /LAB42/G 

COMMON /LA B47 / EPPS1 s EpS2 ,EPS3,EPS4» ITIR 

15 FORMAT (* COUNTER NOT SUFFICIENT FOR MINIMISATION*) 

16 FORMAT (* NORMAL CONVERGENCE AFTER* ,13,* IT IRATION*) 

24 Format <* STEP=*,E12»2,* C0SINE=*,F9°6>* I D I=*,F6®2) 

■25 FORMAT (* D I RECT I ON* , 5 E 1 1 ® 5 , / » 10X , 5E1 1 0 3 ) 

26 FORMAT (IX, *VAR I ABlES* > 10F6 * 3 ) 

28 FORMAT (* MINIMISATION OVER R=*»F11„3) 

29 FORMAT (* GRADIENTS*, 5Ello3,/,10X,5Ell o 5) 

51 FORMAT (* STEP LENGTH LESS THAN OoOIX-MIN WAS OBTAINED FOR N TIMES* 
1*,/,* ST OPING MINIMISATION FOR THIS R AT *,I5,* ITIRATIONS*) 

PRINT 28, R 

PRINT 26, (X( I ) ,I=1,N) 

L7 = 0 
LI = 0 

CALL PENAL(FC»FP»FT »X»N»R»GX) 

IF(FT®GTol0oQE20) GO TO 17 
CALL GRAD ( X,N,G»R) 

PRINT 29, (6(I),I=1»N) 
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DO 5 13 = 1 » N 
DO 5 12 = 1 , N 
H{ 12*13) = 0 o 0 
i H( 13,13) = 1,0 

1 DO 6 II = I , N 
D ( 1 1 ) = 0 c 0 
DO 6 12 = 1 , N 

=> D ( 1 1 ) = D ( 1 1 ) — H ( 1 1 » 1 2 ) #G ( 1 2 ) 

PRINT 25, (D(I),I=1,M) 

SO = 0 o 0 

CXCX = 1 o 00—4 
DO 10 I = l,N 
? n SO = S9+D( I ) *G ( I ) 

IF(S9oLT„0,P) GO TO 2 
DO 23 12 = 1 , N 
IF (H( 12,12) .NE. 1.0) GO TO 1 
23 CONTINUE 

PRINT 132, LI, R 

132 FORMAT ( 15X,*DERIVATIVE OF FUNCTION ALONG DIRECTION OF *,/,15X- 
1*LIKEAR MINIMISATION IS NOT NEGATIVE AT*,IA,* I T I RAT I ON FOR R = *» 
2F 10 <>6 ) 

RETURN 

2 CONTINUE 
XMM = 1.0 

DO 199 I = 1 » N 

IF ( XMMoGToX ( I ) ) XMM = X ( I ) 

199 CONTINUE 

XMM = XMM/100o0 
CALL SEARCH(SIG»X,N*R) 

IF (SIGoLT.XMM) L7=L7+1 
IF ( L7 o GT 0 5 ) GO TO 50 
DO 7 11= 1 »N 
X(I1) = X( I 1 ) +S I G*D (II) 

7 G 1 ( T 1 ) = G ( 1 1 ) 

LI = L 1+1 

IF (Ll.GT.TTIP) GO TO 14 
CALL PENAL (F0»FP,FT1,X*N»R»GX) 

FT = FT1 

CALL GRAD ( X » N »G » R ) 

PRINT 29, (G(I),I=1,N) 

ET = AMAX1 ( 1.0, SIG) 

L 2 = 0 

DO 3 I = 1 , N 

IF(ABS(ET*D(I) J.LT.0.001) GO TO 3 
L2 = L2+1 

3 CONTINUE 

I F ( L2 o EQoO ) GO TO 13 


SI 

=T 

0.0 

S2 

= 

0.0 

S3 

ss 

0 o 0 

DO 

11 

I = 1 , N 

SI 

= 

S 1+G ( I )**2 

S2 

= 

S2+D ( I)**2 

S3 

= 

S3+G ( I 5 *D ( I ) 

S2 

2= 

SORT ( S2 ) 



SI = SORT (SD 
CS = ABSIS3/S2/S1 ) 

• . PRINT 24?SIG?CSjS2 

PRINT 26,(X(I)»I=1»N) 

IF (SIGoLE,10 o 0E-6) GO TO 1 
IF ( Si o LE oCXCX oOR o S2 &LE oCXCX ) GO TO 13 
151 DO 9 II = 1»N 
9 Y ( 1 1 ) = G ( 1 1 )-Gl( 1 1 ) 

SUM2 = OoO 
DO 156 I = ] ,N 
SUM 1 = OoO 
DO 157 J = 1»N 

157 SUM1 = SUM! +X ( J ) *H ( J » I ) 

156 SUM 2 = S'JM2+X< I )*SUM1 
A = SUM2 

IF(AoLToOoO) GO TO 1 
CALL HMATRI (HsSIGsN) 

G" 1 TO S 
PRINT5UL1 
GO TO 17 

14 PRINT 15 

GO TO 17 

13 PRINT 16 » LI 

17 RETURN 

END 


SIBFTC HMAT 

SUBROUTINE HMATRI (H,SIG»N> 

COMMON/LAB23/D 

COMMON/LAB38/Y 

DIMENSION H5(ll)»H9ai»ll).HlO(ll,ll)*HS(ll> 
DIMENSION H6{ll,ll)fH(ll»ll)fY(ll)»D(ll)»H4{lir 

Hl*0o 

■DOIT * 1 » N 
HI = H 1+D ( I )*Y( I ) 

H 5 ( I ) = OoO 
H4( I ) * OoO 

DO 1 J = 1 » N 

H4( I ) = H4( I )+H( I , J >*Y ( J ) 

H5( I) = H5( I )+H( J,U*Y(J) 

]. H9( T*J) = D ( I )*D(J) 

HU = o 

D02I = 1»N 

H 1 1 = Hll+YI I )#H4( I ) 

DO 2 J = 1»N 

2 H 10 ( I »J) * H4( I )*H5( JJ 
DO 3 I = 1»N 

DO 3 J = UN 

H6(I»J) = H(I»J)+SIG/H1*H9(I»J)-H10(I»J)/H11 

3 CONTINUE 

DO 5 I * 1»N 

DO 5 J * 1 » N 

5 H ( I , J) = H6 ( I » J ) 

RETURN 

END 



Z 6 = IZG 

IF ( { ZS( I »J)-ZP) oG'EoOo5)ZP = ZP + l o 0 

I F f ( ZS ( 1 + 1 » J 1 )**ZG ) oGE *0 ® 5 ) ZG = ZG+l o 0 
ZS(I»J) = 7P 

ZS( I+1,J1) = ZG 
TRMS = ZP/ZG 

C ESTIMATION OF TOTAL CORRECTION »XT=2X 

0=AM( I,j)*(Zp+ZG)*C0S(ALFA)/f2»0*CEN<I> ) 

A3=ATAN ( SORT ( l o 0-G**2 ) /Q) 

XT= (TAN( A3 )-A3-TAN( ALFA)+ALFA)*{ ZP+ZG ) / < 2.0*TAN ( ALFA ) ) 
PRINT 8 9 XT 

IF { I COMP ( I » J ) oN E c 0 ) GO TO 43 

C DISTR1 BUT I ON OF TOTAL CORRNo BETWEEN PINION AND GEAR 

CALL CORREC (ZP 9 ZG 9 TRNS 9 XT 9 XP 9 XG) 

CORR ( I 9 J ) = XP 
COP. ( 1+1 9 JI ) = XG 
14 = I COMP ( I + l.Jl) 

IF(I4gMEoO) CORR ( 1 + 1, 14) = CORRII+UJl) 

GO TO 44 

•3 CORR ( I +1 > J 1 ) = XT— CORR ( I 9 J ) 
n '-n'T 47,J,I 

-• J. r 1 :\>T 762 9 ZP 9 XP 9 ZG 9 XG 
C~2 = 1+1 9 JI )*ZG 

DPI = AM(I,J)*ZP 
V = 3 s 1416*SPD< I 9 Ji*DP 1/60000 oO 
XI = C0RR(T 9 J) 

X 2 = CORR ( 1 + 1 ,J1 ) 

PRINT I 3 R 9 I 9 J 
J M = I + 1 

C CALCo OF M'-MEASUREMENT DATA 

CALI. INSPEClZp»DPl*Xl»V*AMd«J) 9 IOLTY) 

PRINT 1399TN9J1 

CALL INSPEC(ZG,DP2.X2,V9AM( I 9 J) 9 IQLTY ) 

41 CONTINUE 
40 CONTINUE 

DO 100 I = I 9 NGRP 

II = IP C I ) 

DO 100 J * 1 9 II 
J1 = IP{I+1)+J 

ICO TR(I.J) = ZS{I*J)/ZS(I+1»J1) 

12 = 0 

13 = 1 

FSPD(l) * RPMI 
FSPDL(l) = FSPD(l) 

DO 120 I = 1 » MGRP 
II = I P C I ) 

14 = 13 
13 = 0 

DO 130 J = 1 9 1 1 
DO 130 K = 1»I4 
13 * 13+1 

FSPDL ( 13) = FSPD(K)*TR(I ,J> 

130 CONTINUE 

IF ( I « NEo TBGSHi GO TO 140 
DO 150 M = 1914 
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15" FSPDP(M) = FSPD(M) 

N1 = 14 

1 4° I F ( ( I-l ) 0 NFo IBGSH) GO TO 160 
DO 170 M = 1,N1 

13 = 13+1 

170 FSPDKI3) = FSPDP(M) 

160 DO 180 M = 1,13 
180 FSPD(M) = FSPDL(M) 

120 CONTINUE 

14 = 13 
201 X = 0 o 0 

I = 1 

218 IF ( X oGEoFSPDL ( I ) ) GO TO 219 
X = FSPDL(I) 

N = I 

219 I * 1 + 1 

IF { I o LEo I 3 ) GO TO 218 
FSPDL<N) = FSPDL(I3) 

F.SPDL ( 1 3 ) = X 
13 = 13-1 

IF ( 1 3 o GT o 1 } GO TO 2Ql 
PRINT 401 
J = 0 

DO 402 I = 1,14 

IF (I cEQo 141 GO TO 403 

IF ( INTtFSPDUI ) ) oEQo INT(FSPDL( I+UU GO TO 402 
4 '~2 J = J+l 

PI = RPML*PHI**< J-l ) 
d = (Pl-FSPDL(n ) /Pl*100 o 0 

PRINT 404, J,P1 ,FSPDL( I ) ,P 
402 CONTINUE 

SI GAL = 15,o 
DO 45 K = l.NGRP 
LI = IP ( K ) 

DO 46 L = 1 , LI 

TORQ = 716„0*HP*ETA**K/SPD(K»Lj*1000.0 
IF ( TORQoGT oTMAX) T0RQ=TMAX 
TLD = TORO/DSH ( < ) /2 o 0 
AREA ( K »L ) = TLD/SIGAL 

AREA(K+1,L1) = AREA(K,L)*DSH(i<+l )/DSH(K) 

46 CONTINUE 
45 CONTINUE 

PRINT 398 , ( ( AREA ( K > L ) »L*1»N6MAX) ,K=1,NSH) 

CALL SHBG 
CALL CLUTCH 
RETURN 
END 


SIBFTC CORREC 

SUBROUTINE C0RREC(2P,ZG»TRNS,XT ,XP,XG) 

DIMENSION XM ( 12 ) ,SLP(12) .XMK13) »SLP1( 13) 

DATA XM/— o ,35 , — 0 o 2 5 ,-0 o 14 »*-0 0 05 » 0 o 08 »0 c 17 ,0 » 26 ,0*38 »0 0 50 ,0 
10.72»0o84/ 

DATA SLP/- ,8935 7535 s- ,6085 »- 0 4835 »- 0 3775 *- 0 2735 ,-o 1762* 

1 — o 0 043 6 , + o 05 67 » + * 1 3 5 1 , + o 1942 / 


62, 

O 0906, 



DATA XMl/-0c38»“0ol8,-*0cl29-0o06,0 o 00»0on7,0ol5,0o23,0o34,0o44, 
I0o57s0o72»0o90 / 

DATA SLP1 / “•0o3060,“0o0559,«0o0^36 ,-0o0l74 5 0 o n,0o0332,n o n663, 

lOo 1050 »0o 1620* Oe 2325,0 o 3305. 0« 4750 *0o 6850 / 

X = XT/2 o 0 

C DISTRIBUTION OF CORRECTION 

ZZ= (ZP+ZG ) /2 o 0 
IFtTRNS oGTo 1 o 0 ) GO TO 731 
C FOR SPEEDING DOWN PAIR OF GEARS 

XX1=SLP( l)*(ZZ-40 o 0)/50.0+XM(l)-X 
DP 735 N=2 » 12 

XX2=SLP(N)*(ZZ-40»0)/50 o n+XM(N)-X 
I F ( XX 1*XX2 ) 736,736,734 

734 X X 1 = X X 2 

735 CONTINUE 

736 SL=SLP(N-1 ) -XXI* ( SLP ( N ) -SLP ( N-l ) )/(XX2-XXl) 

GO TO 751 

C FOR SPEEDING UP PAIR O c GEARS 

731 I F ( X oGEo (— 0 0 3 ) ) GO TO 732 
I F ( 7Z oGTc 5 0 o 0 ) GO TO 739 

732 XX1*SLP1 (1)*(ZZ-4G,0)/50„0+XM1(1>-X 
DO 795 M=2,13 

XX2=SLP1 (N)*(ZZ-40a0)/50 o O+XMl(N>»X 
I c ( XX 1*XX2 ) 737,737,738 

738 X X 1 = X X 2 
795 CONTINUE 

737 SL = SLPl ( N-l )~XXl* ( SLP 1 ( N )— SLP 1 C N-l ) ) / ( XX2-XX1 ) 

751 XP=SL*tZP-ZZ)/50oO+X 

XG*XT**XP 
GO TO 761 

739 XP=X 
XG = X 

761 RETURN 
END 


Tt I B F T C I N SPEC 

SUBROUTINE INSPEC(Z,DP ,X,V,AM,IQLTY> 

C Z = N 0 o OF TEETH, X=CORRECTION 

C CALC, NOo OF TEETH, AN OVER WHICH MEASUREMENT IS TAKEN 

135 FORMAT ( 1 DX ^MEASUREMENT OVER* , 14, 3X ,*TEETH*/ 

1 10X ,*N0MI NAL VALUE=*,F9o3 ,3X,*MM*,/ 

2 10X,*LOWER LIMIT =*,F9»3,3X ,*MM* , / 

3 10X,*UPPER LIMIT=*,F9 0 3,3X,*MM*,/I 

823 FORMAT ( IX , *CAUT I ON— SURRFACE SPEED EXCEEDS 30M/S*) 

I F ( ABS ( X ) oGE« 0 o 4 ) GO TO 810 

AN=Z*20o0/180 o 0+0o5 

GO TO 811 

810 A=Z/3 o 1416*SQRT(0 o 132+4o526*X/Z+4o526*(X/Z)**?)-(Z+2 o 0*X)*0 o 116 

AN = Z*20 o OX 1 80 o 0+0 o 5+a 

811 N=AN 
RN = N 

IFHAN-RN) oGTo 0 o 5 ) GO TO 812 
AN = RN 

GO TO 815 ' 

812 AN=PN+lo 0 

815 ALFA=20o0*3 0 1416/180o0 

C CALCo THE nominal value of measurement 

BT=AM*(Z*C0S (ALFA)* (TAN (ALFA) -ALFA) +<AN-Oo 5 )**, 141 6*COS ( ALFA) 


C 


1+2 oO*X*SIN( ALFA) ) 

SELECTING TOLERANCES DEPENDING ON QLTYo 
I F ( DP a LEo 50 c 0 ) GO TO 101 


AND DIAMETER OF GEAR 
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GO TO 815 
812 AN=RN+loO 

815 ALFA=20oC*3ol416/13Co0 

CAL Co THE NOMINAL VALUE OF MEASUREMENT 

3T*AM* (Z*CDS( ALFA)*! TAN (ALFA) -ALFA )+( AN-0«5 ) *3 0 1416*C0S { ALFA) 

1+2 =~*X*SIN( ALFA) ) 

SELECTING TOLERANCES DEPENDING ON QLTYo and DIAMETER OF GEAR 


IF( DPuLE. 

■ 50 o 

0) 

GO 

TO 

101 

I F ( DP t- LF ; 

,noo 

) 

GO 

TO 

104 

I F ( D P o L E « 

: 2 10o 

0) 

GO 

TO- 

106 

IF( IQLTY 

c ‘ 0 c 

8) 

GO 

TO 

109 

If ( IOLTY 

o EQc 

7) 

GO 

TO 

110 

I F ( IQLTY 

o :0 c 

6) 

GO 

TO 

111 

I F ( TOLTY 

o H.Q O 

5) 

GO 

TO 

112 


109 T0LL=-75oC 
TOLU=-45 c 
IOLTY=8 

GO TO 130 

110 T OLL = -62 o 
T OLU=-41 o 
I 0LTY=7 
GO TO 130 

111 T OLL = — 60 o C 
T 0LU=— 45 o 

I 0 L T Y = 6 
GO TO 130 

112 T 0 L L = - 5 2 o 
T0LU=-41 o 
I 0LTY=5 
GO TO 130 

1~ + IFUOLTY .cIQo 8) GO TO 105 
IF( JOLTY c EO o 7) 60 TO 107 
I F ( JOLTY oEOo 6) GO TO 115 
tF(TQLTY cZQo 5) GO TO 117 
105 T OLL=-85 v 
T 0LU=— 5 1 o 
I OL T Y = 8 
GO TO 130 
107 TOLL=-70o 
T 0LU=— 47 o 
I QLTY=7 
GO TO 130 
115 T OLL--68 e 
TOLU=-51 o 
I QLTY=6 
GO TO 130 
117 TOLL=~5fe o 
T0LU=-49o 
I OLTY=5 
GO TO 130 


106 

I F ( IOLTY 

c EQ o 

8 ) 

GO 

TO 

119 


IF ( IQLTY 

oEQo 

71 

GO 

TO 

121 


IF( IOLTY 

o EQ o 

6) 

GO 

TO 

120 

119 

IF ( IOLTY 
TPLL=-94o 

o EQ o 

51 

GO 

TO 

122 


TOLU=-56o 
I OLTY=8 
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GO TO 130 

121 TOLL=~79o 
T 0LU=— 53 o 
IOLTY=7 
GO TO 130 

120 TOLL=**?5 o 
TOLU=**56„ 

I O L T Y = 6 
GO TO 130 

122 T,0LL = 85o 
T0LU=«68o 
I OL TY = ? 

GO TO 130 

113 = n lNT 823 
130 -■ T L =B T +T 0 L L*0 „ 0 0 1 
3TU=eT+TOL'J*O o O01 
N7 = A.N 

PRINT 135jNZ»BT jBTL>BTU 

RETURN 

END 


SIBFTC SHBG 

SUBROUTINE SHBG 

DIMENSION AM<6»6)»TR(6»6) »CEN (6) 9 AAMOD ( 2 A) 9 SPD (6 »6 ) »SPDM(6»6) 9 
IIP (6) s DSH ( 6 V ? I COMP ( 6 s 6 ) »ZS ( 6 *6 ) » AREA ( 6 »6 ) 9 ALSHC 6 ) »ANG<6 ) 9 DI S { 6 > 6 ) , I CL ( 
21 CL (3 ,6) 9 1 CAR< 6 ) ,CORR<696> 

COMMON/LA EC/RLL 9 PLR »BEN >EQ I VM 9 BDCLHjBDCRH, DEFT, RP 

COMMON /LAB 2 / NSPD 9 NSH , I BG SH 9 I L AP , NGMAX 

COMMON/LAB3/HP 9 ETA 9 RPMI 9 RPML » RPMH , TMAX ,PHI 

COMMON/LAB5/IP»ZS,SPD 

COMMON /LAE 6/ I CL 9 ICAR 

C0MM0N/LAB7/SIGW»FSC.FNUT*FS»E,ALH 

COMMON/LAB A/ PLF.CUTF* DSP 9 DSJ 

C0MM0N/LAE9/AM»TR9DSH,ALSH 

COM M ON/L AE 10 / ANG 9 D I S 

C0MM0N/LAE13/ ICOMP 9 SPDM 

COMMON /LABI 5/ CCRR 

COMMON/LAE AA/CEN 9 AREA 

COMMON /LAB A6 / AAMOD 

175 FORMAT ( 15X s *DlAo OF SHAFT EXCEEDS ALLOWABLE LIMIT*/) 

177 FORMAT (1X,*DATAS FOR SHAFT* 9 1 3 9 / 9 IX 9 *MAX I MUM DIAMETER** »F6» 0 , 

1/ 9 IX 9*MAXI'TJM DEFLECTION**, F6o3»/»1X**RADIAL LOAD LEFT BR=*9Fo o 09/ 

3 9 1 X 9 *P. AD I AL LOAD RIGHT Br=*,F9*0,/1X»*STATIC CAPACITY LEFT BR=*,F9°0 
Au 09 /9lX 9 *STATlC CAPACITY RIGHT BR=* ,F? .0 ,/lX 5 *DYNAMIC CAPACITY LEFT BR 
5T BR=*,F9,0,/9lX9*DYNAMlC CAPACITY RIGHT BR=* , Fo c 0 9 / 9 1 X 9 *AVArAGE D 
6DYNAMIC CA^o LEFT BR=*9F9 oO,/* AVRAGE DYNAMIC CAP.RT BR=*<Fq„09/ 

79 1X,*MAXIM'JM BENDING MOMENT** ,F9 «0 */* IX, *MAXIMUM EQIVALLENT BENDING 
8 G MOMENT =* 9 F9 c 0 ) 

NGRP = NSH -1 

RPM I N = 716 « 0*HP*ETA**NSH/TMAX 
C DESIGN STARTS FROM SHAFT 1 

ALFA = 3« 1416/20 oO 
L = 1 

101 RLMXL = 0 o 0 
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RLMXR=0« 0 
EQMAX-0 o 0 

BM AX = 0 „ o 
DCLM= 0 <j n 
0C D M=0 o 0 
DEFL = OoC 
A D C o = 0 o 0 
ADCL = "*oO 
9 AWC=0 0 o 

I F ( L « N E o 1 o AN D o L j N E o NSH ) GO TO 10 
IF(LoNEol) GO TO 11 
I 1 = 1 

I 2= I P ( 1 ) 

D I S1=DSP 
F1=PLF 

IF(ANG(1) oGTolSOoO) GO TO 13 
AM = ANG ( 1 } *1 10 o 0 

G 0 To 12 

13 AM = ANGddUOoO 
GO TO 12 
3 1 11 = 1 

I2=IP(NSH*1| 

D I S 1 = DSJ 
FI = CUT F 
A N = 0 c 0 

12 DO 14 1=11,12 
AMI = AM ( L ♦ I ) 

RP = SPD(Ld ) 

I F ( SPD ( L 9 1 ) o LT RPMIN) RP=RPMIN 
TPPQ=716c*HP*ETA#*L/SPD (L»I ) *1000,0 
IF (TOROoGT 0 TMAX> TORO=TMAX 
F2 = 2oO*TOPO/ZS(L»I)/AM(Ld )/0o937 
XL=ALSH ( L) 

DIS2=DIS(L»I ) 

CALL SINGLE < FI , F2 , XL , D I S 1 » D I S2 » AN ,T0R0 » ALH > 
IF(RLLoGT 0 PLMXL) RLMXL=RLL 
IF(RLRcGToRLMXR) PLMXR=RLR 
IF(FQIVMoGToEOMAX) FQMAX=EQIVM 
IF(BENoGTo3MAX! BMAX=BEN 
I F ( BDCLH o GT o DCLM. ) DCLM=BDCLH 
IF(BDCRHcGToDCRM) dcrm=bdcrh 
IF(DEFT oGToDEFL) DEFl=DEFT 
ADCR=ADCR+BDCRH 
ADCL=ADCL+BDCLH 
ANC=ANC+lcd 
14 CONTINUE 
505 STAL= RLNXL*2oC 
ST AR = R LMXR *2 o 0 

C CALCo SHAFT DIA. FOR MAXMo EQVo BoM.oAND CHECKING DEFLECTION 

D I A* ( 32oO*EOMAX*FS/(3.1416*SIGW*FSC/FNUTJ >**0,3333 
173 DEF = DEFL/ ( E*3ol416*DIA**4/64 o 0) 

I F ( DEF oLT, 0,. 01 5* AMI ! GO TO 171 
DIA = DIA+loO 

I F ( D I A o GT e 1 00 e 0 ) GO TO 174 
GO TC 173 

DIA1 = (EO'1AX/2oO)**Oo3333 


171 
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IF(DIAlvGT.DlA) DIA = DIA1 
■“DIA = D T A 
Rt l ~ ,M H I A 

■'SH (L ) =DI * 

T F(|_ ,NE, MSH) GO TO 176 
DI -=1 ,06*01 i 
DSH(L) = DIA 
GO TO 176 
174 POINT 175 
176 ADCR=ADCR/ANC 

adcl=adcl/anc 

PRINT 177 > L 9 DI A > DEF >R.LMXL > RLMXR > STAL 9 STAR » DCLM »DCRM * ADCL » ADCR » BMAX 
1 s RQMAX 
L = L+1 

IF(L =LEo MSH) GO TO 101 
GOT 0 99 

C CAL CM o STARTS FOR SHAFTS OTHER THAN 1ST, AND LAST ONES 

10 N 1 « I P ( L ) +1 

N2 = I P ( L ) +1 P ( L~1 ) 

12 = I P < L ) 

DO 16 J * N 1 » N 2 
A w 1 = AM ( L 9 J ) 

AN = ANG(U 
RP = SPDCLsJ) 

IF (RF 0 LT C RPMTN) RP = RPMIN 
TORO = 716 o D*HP*ETA**L/SPD(L»J)*1000o0 
I F ( TOROo GT o TMAX ) TORO=TM AX 
FORC = 2oC*T0PQ/ZS<L»J)/AV<L*J)/0e937 
C CALC » FOR ^ARTICULAR SHAFT FOR VARIOUS LOADINGS 

DO 16 I = 1*12 

FORCE = 2oO*TORQ/ZSIL»I)/AM(L9l )/0o937 
I F ( AN— 180 tO ) I 79 I 89 I 8 
18 AN = AN+4C,0 
GO TO 15 

17 AN = AN-40,0 
15 XL = ALSH(L) 

C CHECKING I c ANY FORCE IS OUTSIDE THE SPAN BETNo THE BEARINGS 

I F ( DIS (L » I ) oLToXLoANDoDIS(L» I ) » GT 0 0 o 0 * AND „ D I S ( L » J ) °LT 0 XL 0 AND 0 DI S ( L 9 J ) c 
1 » J ) oGT.O.C) GO TO 19 
I F ( DIS(L9 I ) cLToOoO) GO To 154 
IF(DIS(L»I ) 0 GT 0 XL) GO TO 155 
IF(DIS(L»J) oLToOoO) GO TO 156 
IF(DIS(L»J) 0 GT 0 XL) GO TO 157 

154 DIS1*DIS1L»I > 

DIS2*DIS(L»J) 

F1=F0RCE 

F 2 = F 0 R C 
GO TO 159 

155 DIS 1 =DIS(L*I)-XL 
DIS2=DIS(L»JJ 

p 1 = FOP CF 
F2=F0RC 
GO TO 159 

156 DIS1=DIS(L »J) 

DIS2 = DIS(L9I > 
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Fl-FORC 
F 2 = F 0 R C E 
TO 159 

157 DIS1=DIS(L*J)-XL 
DIS2 = DIS(L , I ) 

F 1 = FOR C 
F 2= FORCE 

159 CALL SINGLE ( FI , F2 , XL , D I Si » D I S2 , AN , TORO , ALH ) 

GO TO 168 

C CALC* FOR GEAR FORCES ALL IN BETWEEN THE BEA9INGS 

19 FOP CA = FORC*COS ( AN/ 5 7o 3 ) 

FOR CP = FORCES I N < AN/5 7 0 3 J 

RRl=-( 01 S (L» I)*FORCE+DlS{L, J) *FORCA) /XL 
RLl=-( FORCE+FORCA+Rf?l ) 

RR2=**FORCP*D I S ( L ? J ) /XL 
RL2=— < FOPCP+RR2 ) 

RLR=SQRT l PR1**2+RR2**2 } 

RLL=SQRT< RL1**2+RL2**2 ) 

IF(DIS(L»n-DlS(L*J) ) 21,22,22 
21 '<1 = 1 

F\/A = FORCA 
FVP=FORCP 
ARO=DIS(L,H 

Dl=(DIS(L,.J)-DIS(L»m/10 o 0 
CO TO 23 

22 <L=0 

FV A = FO r ?CE 
F V D = o 0 0 
A P '!=D I S ( L 9 J) 

Dl= ( D I S ( L , I >*DIS(L,J> )/10o0 

23 BEN = C c 0 

C CALCo BoMoS AT 10 SECTIONS OF SHAFT AND SELECTING THE MAXMo B.Mo 

DO 24 M=1 » 10 
A = M 

A R M = A R M*f D 1 

BENX=SQRT ( (RL1*ARM+FVA*A*D1 >**2+(RL2*ARM+FVP*A*Dl >**2 > 
IF(BENX-BEN) 24,25,25 

25 BFN=BENX 

24 CONTINUE 

EQI VM=SQRT ( (TORQ/2oO)**2+BEN**2> 

BDCRH = {ALH*60 o 0*RP)**0 o 3333*RLR/100 o 0 
BDCLH = (ALH*60o0*RP)**0o3333*RLL/100o0 
C CALCo DEFLECTION AT GEAR POSITIONS 

XLL=XL**2 
DEN0=6 o 0*XL 
B1=XL«DIS(L, I > 

B2=XL“DI S ( L , J ) 

D 1 1 = F0RCE*D I S ( L » I )*B1* ( XLL-DIS ( L , I )** 2-8 1**2 ) /DENO 
D 2 2 = F0RC'*D I S ( L , J )*B2# ( XLL-DI S ( L » J} **2— B2*-*2 ) /DENO 
IF ( DIS (L , I ) oLEoDIS ( L, J ) ) GO TO 26 

D21=FORCE*DlS ( L , J ) *B1* ( XLL-DI S ( L , J ) **2~B 1**2 ) /DENO 
D12 = D2V* FOR C/FORCE 
GO TO 27 

26 Dl2 = FORCED I S ( L » I ) *B2* ( XLL-DIS ( L , I ) **2-B2**2 ) /DENO 
D21 = D12*FORCE/FORC 

27 Dl= ( D1 l**2+Dl2**2+2oO*Dll*D12*COS( AN/57 *3) ) **0e5 
D2= ( D2 2**2+D2l**2+2o 0*D22*D21*COS( AN/57 *3 ) ) * # C*o 5 
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DEFT = AMAX 1 ( D1 » D2 ) 

150 CONTINUE 

C CALC, MAXMo DEFLNo BM«,EQVo BMoFOR SHAFT DYN® CAPACITY AND RADIAL 

C LOADS FOR BEARINGS 

168 I F ( DEFToGT , DEFL ) DEFL=DEFT 
IF(BDCRHoGToDCRM) DCRM-BDCRH 
IFCBDCLHoGToDGLM) dclm=bdclh 
I F C BFNoGT,3MAX) BMAX=BEN 
TF( EfJlVMoGTu EQMAX ) EQMAX=EQI VM 
IF(RLRoGToRLMXR) RLMXR=RLR 
IF(RLLoGToRLMXL) RLMXL=RLL 
ADCP=ADCR+3DCRH 
ADCL=ADCL+BDCLH 
ANC=ANC+lo3 
: j CONTINUE 
CO TO 505 
99 RETURN 
END 


SIBFTC SINGLE 

SUBROUTINE SINGLE (Fl»F2»XL»DISl»DIS2»AN»TORQ»ALH) 

common/labc/rll*rlr»ben*eoivm*bdclh »bdcrh»deft ,rp 

C CALCc WHEN ONE LOAD IS OUTSIDE THE SPAN B E T N • BEARINGS 

B 1 = XL-DIS2 
BEN1=F1*DIS1 
RLL2 =-F 2*3 1 /XL 
RLR2=~( F2+RLL2) 

BEN2 = RLL2*DIS2 
T 1 =-Fl*(A3S(DISl )+XL) /XL 
T 2 = F1*ABS( DIS1WXL 
IF(DISloLToOoO) GO TO 200 

RL.L * ( T2**2 + RLL2**2+2 o 0*T2*RLL2*C0S(AN/57 o 3) ***0 o 5 

PLR = ( Tl**2+RLR2**2+2 o 0*Tl*RLR2*C0S{ AN/57 0 3 >^*0 o 5 
CO TO 201 

200 PLL *(T1**2+RLL2#*2+2 o 0*T1*RLL2*C0S(AN/57o3) ***0,5 

RLR =(T2**2+RLR2**2+2„0*T2*RLR2*C0S(AN>57 o 3J )**0 o 5 

201 BEN = ( BEN 1**2+ BEN2**2+2 0 0*BENl*BEN2*COS { AN/57 e 3 ) )**0,5 

EQIVM = SORT ( (T0RQ/2 o 0) **2+BEN**2 ) 

DFFl = F2*31*DIS2*(XL**2“DIS2**2— Bl**2)/6 o 0/XL 
DEF2 = BENl*DlS2*Bl/4.0 

DEFT =(DEFl**2+DEF2«*2+2 o 0*DEFl*DEF2*CCS(AN/57.3) '**0®5 

BDCRH = ( ALH*60 o 0*RP)**0 o 3333*RLR/100 o 0 

BDCLH = ( ALH*60 o 0*RP ) **0 o 3333*RLL/100 o 0 

RETURN 

END 


$ I BFTC CLUTCH 

SUBROUTINE CLUTCH 

logical nocl 

COMMON /LAB2/NSPD »NSH? I BGSH > I LAP > NGMAX 

C0MM0N/LAB3/HP,ETA»RPMI »RPML»RPMH,TMAX,PHI 
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C Q M M 0 M / L A 5 4 / G D C H , 6 D P , R H 0 
COMr-*ON/LAB5/IP»ZS,SPD 
C0MM0M/LAB6/ICL* ICAR 
C0MM0M/LAB9/AMfTR»DSH»ALSH 
COMMON /LAE 10/ ANG.DIS 
C0MM0N/LAB12/N0CL 
C0MM0N/LAE13/IC0MP» SPDM 
C0MM0N/LAB15 / CORR 
C0MM0M/LAB44 /CEN , AREA 
COMMON /LAB46/AA MOD 

DIMENSION AM(6?6)»TR(6»6) ,CEN<6> »AAMOD(24) ,SPD(6»6) *SpDM(6*6) » 

IIP (6) , DSH ( 6 ) ? I COMP (6, 6) »ZS( 6,6) , AREA (6,6) »ALSH(6) ,ANG(6) *DIS(6*6)»ICL<3 
2 I CL (3, 6) ? I CAR ( 6 ) ,C0RR(6,6) 

DIMENSION GDSH ( 6 ) ,GDG(6,6> ,GDA(3,6) »GDB ( 3 ,6 ) » AO ( 3 ,6 ) , GD ( 6 ) » FW { 6 *6 ) 
DIMENSION TRQ(3 ,6) ,AOI (4) »GDAI (4> ,GDBI (4) , M ( 3 , 6 ) , EQ I V ( 3 , 6 ) ,RPM { 6 ) . 

IF (NOCL) RETURN 

DAT A AO I j GCA I , GDB I / 950 ° 0 » 1440 . 0 » 2200 o 0 , 4200 o 0 ,0 o 002 13 , 0 » 0 06 >0 ° 02 13 
1 » 0.046 »0c Cl 125 *0*0256. 0.0627 *0*14447 / 

ASSUMED CLUTCHES EXIST ONLY ON SHAFTS 1 AND 3 
105 F0RMAT(//1X.*DESIGN DETAILS OF CLUTCHES*/) 

179 FORMAT (1XjI2,4X9F6o3»4X 3 F6o394X,F9 »2) 

45“ FORMAT (IX,* SWITCHING POWER PER OPERAT ION* ,F3 • 3 ** M— KG* ) 

541 FORMAT ( 1 X »*CLUTCH* , 13 »* IS OF THE T YPE* , 1 3 * / * 1 X ,*AND ITS THERMAL 
1 UTILISATION I S* » F9 o 2 ) 

r - Tc FORMAT (1X»*SHAFT GD2— VALUE INITIAL W1 FINAL W2*) 

716 FORMAT (1X,*DATAS FOR *913.* CLUTCH*.//) 

PRINT 105 
ICT = 0 

TMAX=TMAX/1000.0 
NGRP = NSH— 1 
DO R99 I = 1,3 
DO 999 J = l.NGMAX 
GDA(I.J) =3o 
MU *J) = C 
AOC I , J) = OoO 
999 GDB ( I , J ) = OoO 

C CALCo GD2 OF ALL SHAFTS AND GEARS AND SUMMING THEM UP FOR SHAFTS OTHER 

C THAN 1,2,3 

DO 106 I * 1 ,NSH 

GDSH ( I ) = RHO*3o 1416/8oO*DSH( I ) **4*ALSH < I ) / 10, 0**6 
SUM = GDSH ( I > 

DO 100 J = l.NGMAX 
FW( I , J) = 9 o 0*AM ( I , J > 

GDG( I »J) = RH0*3o 1416/8. 0*FW( I * J ) * ( { AM ( I * J ) *ZS { I * J ) ) **4*DSH ( I ) **4 ) 

1/10 o0**6 

I F ( Io EQ.1.0R.I.EQ*3) GO TO ICO 
SUM = SUM+5DG { I , J ) 

GD ( I ) = SUM 
100 CONTINUE 

106 CONTINUE _ . _ . mlllC . Tr( . 

C CALCo MAX. TORQUES FOR GEARS TO WHICH GEARS ARE CONNECTED 

DO 200 I = 1,3 

DO 200 J = l.NGMAX 

IF(ICL(I»J)oEQoO) go to 200 

TRQ ( I , J ) = 7l6cO*HP*ETA**I/SPD( I ,J) 
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IF(TMAXoLT a TRQ(I > J ) 5 GO TO 305 
TOP = TRQ(I,J) 

GO TO 302 
305 TOP = TMAX 
3 02 CONTINUE 

selecting THE clutches FROM TORQUE consideration 


TF(T0PoLE„10o0> 

GO 

TO 

201 

IF (T0PoLE«.20o0) 

GO 

TO 

202 

IF(TOPoLE.45oO) 

GO 

TO 

203 

IF ( T0PoLEo30o0) 

GO 

TO 

204 


201 Ml I ,J)=1 
GO TO 205 

202 M ( I » J ) =2 
GO TO 205 

203 M ( I » J ) =3 
GO TO 205 

204 MCI » J)»4 

205 MM=M( I , J > 

AO ( I ,J)=AOI (MM) 

GDA( I , J)=GDAI (MM) 

GDB ( I 9 J) =GDBI (MM) 

200 CONTINUE 
GO TO 113 
5 •"'3 MM = M( 14, J4) 

AO ( 14, J4)=A0I (MM) 

GDA( I4,J4)=GDAI (MM) 

GDB ( I4,J4)=GDBI(MM) 

C CHECKING THERMAL UTILISATION FACTORS BEGINS 

113 DO 50 1=1,2 
J 2 = I P ( 1 ) 

K 2 = I P ( 2 ) 

GO TO ( 55,56) , I 

55 J1 = 1 

K 1 = I P ( 2 ) 

GO TO 57 

56 J1 = IPtl) 

K 1 = 1 

C LOOP EXECUTED ONCE ONLY FOR 2ND„ GROUP CLUTCHES 

57 DO 50 J = J1»J2 

IF (IoEQ.l) ICT = ICT+1 

11 = 0 

C MAXMc RPM OF SHAFT 1 

RPM ( 1 ) =RPMI 

I F ( ICL ( 1 » J ) «NEoC ) 11=1 

12 = I P ( 2 ) +J 

I F < ICL(2»I2) oNEoO) 11*2 
I F ( IloEQoC) GO TO 50 
GO TO (58,59) ,11 

58 Ll= 1 

L2 = IP(1) 

C CALCo EQUIVALENTS FOR IDLY RUNNING GEAPS AND CLUTCH PARTS 

DO 60 L = L1»L2 

60 EQ I V ( 1 , L ) = (TR(1,J)/TR(1,L) >**2 

544 FORMAT (2X,*EQUIVALLENT*,51I2*2X*F7.7)) 

GO TO 71 

S' 5 LI = I P ( 2 ) +1 
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L2*IP( 1)+IP(2) 

DO 61 L = L1.L2 

61 EQ I V ( 2 , L ) = (1.0/TR(2,L)/TR(1,J))**2 

71 GD( II) = GDSH(Il) 

DO 72 L * ll»L2 
IF(ICL(I1 »L)oEQo2) GO TO 80 

GD (II) * GD( I1)+(6DA( Il»L)+GDG{ Il,L) )*EQIV( Il,L)+GDB{ I1,L) 
GO TO 72 

80 GD (II) = GD( II )+(GDB( I1,L)+GDG( I1,L) )*EQIV(I1,L)+GDA( I1»L) 

72 CONTINUE 

C MAXo RPM CP SHAFT 2 

RPM(2) * FRM( 1 ) *TR ( 1 , J) 

C LOOP EXECUTED ONCE ONLY FOR 1ST* GROUP CLUTCHES 

DO 50 K = K1 , K2 
IF ( I «, EQo 2 ) ICT = ICT+1 
PRINT 716 » ICT 
13 = 0 

IF( I CL ( 2 » K ) oNEoO) 13=1 

12 = I P ( 3 ) +K 

IFC I CL (3,12) oNEoO) 13 = 2 
IF( I3 .EQoO) GO TO 50 
GO TO (74,75) ,13 

74 L 1=1 

L 2 = I P ( 2 ) 

13 = 2 

DO 76 L * L1»L2 

76 EQI V ( 2 »L ) = (TR(2,K)/TR(2, L))**2 
GO TO 77 

75 LI = I P ( 3 ) +1 
L2* IP ( 2 )+IP ( 3 ) 

DO 78 L = L1»L2 

78 E0IV(3,L)«1,0/(TRC3,L)*TR{2>K) )**2 

13 = 3 

77 GD( 13 )=GDSH( 13) 

DO 79 L = L1,L2 
I F ( ICL( 13, L) oEQ 0 2) GO TO 81 

GD ( I 3 ) = GD( I3)+(GDA( 13, L>+GDG( 13, L) >*EOlVC 13, L)+GDB( 13, L> 
GO TO 79 

81 GD (13) = GD( I3) + (GDB< I3,L)+GDG( 13, L) > *EQ I V { 13 , L ) +GDA { 1 3 , L ) 
79 CONTINUE 
C MAXMo RPM OF SHAFT 3 

RPM ( 3 ) = RPM(2)*TR(2,K) 

C MAXMo RPMS /F SUBSEQUENT SHAFTS 

L = 4 

90 INDEX * ICAR(L-1}*( IP(L-l)-l) 

RPM(L) = RPM(L-1)*TR(L»1*1)*PHI**INDEX 
IF(LoEOoIBGSH) GO TO 91 
92 L = L+l 

IF( LoGToNSH) GO TO 95 
GO TO 90 

91 R PM ( L + l ) = RPM(L-l) 

L = L+l 

GO TO 92 
95 CONTINUE 

DJ0B=10C0o0*100o0/3ol416/RPM(NSH> 
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GD JOB = 3„1416*RH0*5*0*DJ0B**5/8o0/10 o 0*#6 
GD(NSH) = GD ( NSH ) +GD JOB+GDCH 
GD ( 1 ) = GD ( 1 ) +GDP 
5UMGD=0 o 0 

C CALCo OF GD2 OF SHAFT ASSEMBLIES AND SWITCHING POWER PER OPERATION 

DO 150 KK = 1 » NSH 

150 SUMGD * SJMGD+GD( KK)*RPM(KK)**2/7160.0 
PRINT 546 
SINI = 0 * C 
DO 717 IT * 1 »NSH 
PRINT 179* IT»GD( IT) ,SINI *RPM( IT) 

717 CONTINUE 

PRINT 457 * SUMGD 
I F { IcNEol) GO TO 550 
14=11 
J4 = J 

GO TO 551 

550 14=13 

J4 = K+I P ( 3 } 

C CALCo AND CHECKING W 

551 W = SUMGDMO( I4,J4)*100o 0 
PRINT541 ,ICT*M( I4,J4) *W 
IF ( WoGTo 1C Do 0 1 GO TO 209 
GO TO 50 

209 M( I4» J4)=M( I4*J4)+1 

PR I NT 541 , ICT»M( 14, J4) ,W 
GO TO 503 
50 CONTINUE 

TMAX=TMAX*1000 o 0 

RETURN 

END 


SENTRY 



